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ABSTRACT 


We consider the computation of the adjacency characteristic polynomial of a uniform 
hypergraph. Computing the aforementioned polynomial is intractable, in general; 
however, we present two mechanisms for computing partial information about the 
spectrum of a hypergraph as well as a methodology (and in particular, an algo- 
rithm) for combining this information to determine complete information about said 
spectrum. The first mechanism is a generalization of the Harary-Sachs Theorem for 
hypergraphs which yields an explicit formula for each coefficient of the characteristic 
polynomial of a hypergraph as a weighted sum over a special family of its subgraphs. 
The second is a mechanism to obtain the eigenvalues of a hypergraph in terms of 
certain induced subgraphs. We further provide an efficient and numerically stable 
algorithm which combines this information, the set-spectrum of a hypergraph and 
its leading coefficients, to determine the multi-set spectrum of the hypergraph. We 
demonstrate our algorithm by computing the characteristic polynomial of numerous 


hypergraphs which could not be computed using previous methods. 
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CHAPTER 1 


INTRODUCTION 


Spectral graph theory is a well-studied topic which seeks to relate the structure of 
a graph with its spectrum, and vice versa. Results from spectral graph theory have 
found application in combinatorics, computer science, and the social sciences. Most 
notably, in questions regarding the analysis of networks (a famous example being 
Google’s PageRank algorithm). In recent years, there has been increasing interest 
from both academics and practitioners to analyze non-binary interactions between 
agents in a network. Such multi-relational data can be modeled with a hypergraph 
using higher dimensional edges. It is not obvious a priori what definition of the 
spectrum of a hypergraph correctly analogizes the spectrum of a graph. There have 
been various approaches to establish such a definition and we summarize several of 


these approaches below. 


1. Matrices and linear spectral theory. A natural approach is to encode 
the incidence structure of a hypergraph into a matrix and relate the eigenvalues 
thereof back to the structure of the hypergraph. This approach has been widely 
used in applied sciences (e.g., social networks, reaction and metabolic networks, 
protein complex networks, etc. [16]). A benefit of this approach is that one 
can immediately apply linear algebra to study hypergraphs, but considerable 
loss of information about the structure may be incurred when shoehorning the 


hypergraph into a two-dimensional array. 


2. Hypermatrices and low-rank decompositions. One can define a hyper- 
matrix to be a higher dimensional array which maintains desirable properties of 
matrices. In particular, one might seek to analogize the property that the array 
can be decomposed into rank one hypermatrices (i.e., outer products of vectors 
with themselves) as in the classical Spectral Theorem from linear algebra. Im- 
portant contributions to this approach include a generalization of the expander 
mixing lemma [9] and the considerable literature on rank-1 decomposition of 


tensors (e.g., [18]) 


3. Hypermatrices and polynomial maps. Similar to the previous approach 
one considers the entries of a hypermatrix as coefficients of multilinear forms. 
In the case when the underlying hypermatrix is symmetric (i.e., the entries 
are invariant under permutation of the indices) these multilinear forms define 
homogeneous polynomials. The study of symmetric tensors has led to a use- 
ful definition of the spectrum of a hypergraph as the roots of a characteristic 
polynomial which is the resultant of a certain polynomial system whose coeffi- 
cients are drawn from the adjacency hypermatrix, as described in the following 


section. 


We address spectral hypergraph theory from the perspective of hypermatrices 
and polynomial maps. While this approach is admittedly more complex than its two 
predecessors, its merit is that it retains all information of the underlying system. The 
cost of this faithful analysis is paid for in computational complexity: the characteristic 
polynomial of a hypergraph is the resultant of a system of homogeneous multi-linear 
equations. The resultant is an ubiquitous tool in mathematics and has been a concept 
of intense study. Despite this, computing the resultant (and thus, the characteristic 
polynomial of a hypergraph) is a NP-hard in general. In this dissertation we provide 


a new approach for computing the characteristic polynomial of a hypergraph as well 


as a new paradigm for analyzing multi-dimensional networks. 

We begin by providing an analogue of the Harary-Sachs Theorem which yields an 
explicit combinatorial formula for the codegree-d coefficient of a uniform hypergraph 
which is a function of d and the hypergraph’s structure. This formula allows us 
to compute partial information about the characteristic polynomial of a hypergraph 
which is novel as computing the entire polynomial can be intractable. This motivates 
a numerically stable algorithm which allows us to determine the multiplicity of each 
eigenvalue given the set spectrum of a hypergraph and an appropriate number of 
leading coefficients (i.e., an analogue of the Newton Identities). Finally, we prove a 
remarkable property of the spectrum of a hypergraph: the eigenvalues of a hypergraph 
can be written as a union of eigenvalues arising from its induced subgraphs. This 
result is unique to spectral hypergraph theory as it is not true for graphs (c.f. Cauchy 
Interlacing Theorem). We can make further use of this result by appealing to the 
Lu-Man Method to determine the set spectrum of a hypergrpah. It is our hope that 
this dissertation will provide researchers with new tools for analyzing the set and 
multi-set spectrum of a hypergraph. 

The paper is arranged as follows. In Chapter 1 we present the necessary back- 
ground concerning the normalized adjacency spectrum of a uniform hypergraph. We 
then prove an analogue of the Harary-Sachs Theorem in Chapter 2. Using this theo- 
rem we provide an explicit combinatorial formula for low-codegree coefficients of the 
charactersitic polynomial of a hypergraph in Chapter 3. In Chapter 4 we consider the 
computation of the characteristic polynomial given its low-codegree coefficients. In 
doing so we motivate the computation of the set spectrum of a hypergraph which we 
present in Chapter 5. In Chapter 6 we provide several examples of this computation 
as well as examples which highlight the current limits of this approach. We leave the 


reader with open problems in Chapter 7. 


CHAPTER 2 


PRELIMINARIES 


Here we present requisite background maintaining much of the notation of [12]. A 
(cubical) hypermatriz A over a set S of dimension n and order k is a collection of 
n® elements disiz..i;, € S where i; € [n]. A hypermatrix is said to be symmetric if 
entries with identical multisets of indices are the same. That is, A is symmetric 


Signe ee ee a Ee ae for all permutations o of [kK]. An order k dimension n 


symmetric hypermatrix A uniquely defines a homogeneous degree k polynomial in n 


variables (a.k.a. a “k-form”) by 


n 


Fa(x) = S- Qizig.i, Vig Vig +--+ Viz. 


41 29j.05th=1 


If we write x®" for the order r dimension n hypermatrix with 7,,722,...,7, entry 


X1,Uj,...%;, and x" for the vector with 7-th entry x? then the above expression can 


Tr 


be written as 


Ax@k-} = \xh-1 


where the multiplication denoted by concatenation is tensor contraction. Call A € C 
an eigenvalue of A if there is a non-zero vector x € C”, which we call an eigenvector, 


satisfying 
n 
=, k-1 
» Ajig..ip Viz Lig + ++ Lig = AX; . 
12913 ,2:j;2R=1 


Next we offer an important result from commutative algebra to proceed the definition 


of the adjacency characteristic polynomial of a hypergraph. 


Theorem 2.1. (The Resultant, [20/) Fix degrees d,,d2,...,dy. For i € |n|, consider 


all monomials «* (where a is itself a vector) of total degree dj in 11,...,%p. For 


each such monomial, define a variable uj. Then there is a unique polynomial RES € 


Z|{uia}| with the following three properties: 


1. If Fy,..., Fr € Clai,...,¢n| are homogeneous polynomials of degrees dy,..., dn 
respectively, then the polynomials have a non-trivial common root in C” ex- 
actly when RES(F),...,F),) = 0. Here, RES(F,,...,F),) is interpreted to mean 


substituting the coefficient of «* in F; for the variable u;q in RES. 
2. RES(z#,...,0%) = 1. 
3. RES is irreducible, even in C[{uj,.}]. 


Moreover, for i € [n], RES is homogeneous in the variable {uj} with degree 


I @. 


JEln] ,J#t 
Definition 2.2. ([30|]) The symmetric hyperdeterminant of A, denoted det(A) is 


the resultant of the polynomials which comprise the coordinates of Ar®*~!. Let 


be an indeterminate. The characteristic polynomial ¢4(A) of a hypermatrix A is 


pa(A) = det(AZ — A). 


We consider the normalized adjacency matrix of a k-uniform hypergraph, H = 
(V, &). We refer to such hypergraphs as k-graphs and we reserve the language of 
graph for the case of k = 2. For a k-graph H = ([n], E) we denote the (normalized) 
adjacency hypermatriz Ay to be the order k dimension n hypermatrix with entries 

ps ee Lit y%9,.4.5%h) € BCA) 
eens (k-D)! 0 : otherwise. 


For simplicity, we denote ¢(H) = ¢,,,(A) and write 
£ 
o(H) = dood 
i=0 
where t = n(k —1)"~! by Theorem 2.1. Throughout, we make use of the notation 


éa(H) = cq for the codegree-d coefficient of ¢(H). We refer to o(H) = {r : dx(r) = 


0} C Cas the spectrum of H and each A € o(H) is an eigenvalue of H. 


CHAPTER 3 


THE HARARY-SACHS THEOREM FOR HYPERGRAPHS 


An early, seminal result in spectral graph theory of Harary [23] (and later, more 
explicitly, Sachs [33]) showed how to express the coefficients of a graph’s characteristic 
polynomial as a certain weighted sum of the counts of various subgraphs of G (a 


thorough treatment of the subject is given in [3], Chapter 7). 


Theorem 3.1. (/23/,/33]) Let G be a labeled simple graph on n vertices. If H; denotes 
the collection of i-vertex graphs whose components are edges or cycles, and c; denotes 
the coefficient of X”~* in the characteristic polynomial of G, then 
c= XS (1) 20 er CG] 
HEH; 
where c(H) is the number of components of H, z(H) is the number of components 
which are cycles, and |#H C G] denotes the number of (labeled) subgraphs of G which 


are isomorphic to H. 


The goal of this chapter is to provide an analogous result for the characteristic 
polynomial of a hypergraph. The full result is given in Theorem 3.27, but to state it 
here simply: fix k > 2 and let Hy denote the set of k-valent (i.e., k divides the degree 
of each vertex) k-uniform multi-hypergraphs on d edges. For a k-uniform hyper- 


deg =a) 


graph H the codegree-d coefficient (i.e., the coefficient of x of the characteristic 


polynomial of the n-vertex hypergraph 1 can be written 


ca= 2 (-(k-1)") Cu (#H CH) 


HEH, 


where c(/#7) is the number of components of H, Cy is a constant depending only on 
H, and (#H C H) is the number of times H occurs (in a certain sense that is a 
minor generalization of the subgraph relation) in H. 

The quantity (#H C H) is straightforward to compute. However, computing the 
associated coefficient of H, Cy, is more complicated. This notion of an associated 
coefficient of a hypergraph first appeared in [12], where Cooper and Dutle provide a 
combinatorial description of the codegree k and codegree k + 1 coefficient, denoted 
cr and Cr41 respectively, for the normalized adjacency characteristic polynomial of a 


k-uniform hypergraph. 
Theorem 3.2. /12/ Let H be a k-uniform hypergraph. Then 
Ce = —KE(k — 1)" *|E(H)| 


and 


Choi = —C,(k — 1)"-*(# of simplices in H), 
where Cy is some constant depending on k. 


This idea was further studied by Shao, Qi, and Hu where the authors prove 

(restating Theorem 4.1 of [35]), 
ca=(k-1)"™ DP fol€(D)| 
DED 

where D is a certain large family of digraphs, fp is a function of D and €(D) is 
the set of Euler circuits in D. The authors then use their formula to provide a 
general description of Tr2(7’) and Trs3(Z’) for a general tensor T. Our first few results 
are similar to that of [35] (as described in more detail herein), and we use them to 
provide an explicit combinatorial description of Hp and the resulting Cy which yields 


a Harary-Sachs type formula for hypergraphs which is amenable to computation. 


3.1 BACKGROUND 


Our approach relies on the following trace formula for the hyperdeterminant of a 
tensor. In [28], Morozov and Shakirov give a formula for calculating det(Z — A) 
using Schur polynomials in the generalized traces of the order k, dimension n hy- 
permatrix A. Let f : C” — C” be a linear map and let J be the unity map, 


I = (%1,20,...,2n)’ — (21, %0,...,%n)?. Famously, 


log det] — f) = trlogU — f) = 


The characteristic polynomial is defined as the resultant of a certain system of equa- 
tions by definition, so calculating the characteristic polynomial requires computation 
of the resultant. Moroz and Shakirov give a formula for calculating det(Z —.A) using 
Schur polynomials in the generalized traces of the order k, dimension n hypermatrix 


A. 


Definition 3.3. Define the d-th Schur polynomial Py € Z[ti,...,ta] by Po = 1 and, 
for d > 0, 


Pili, So ee 


m=1 d1+---+dm=d 


More compactly, one may define P, by 


exp & tc!) = y Pes 0x. ,tq)z4 
d=1 d=1 


Let f; denote the ith coordinate of Ax®*~!. Define A to be an auxiliary n x n 
matrix with distinct variables A;; as entries. For each 7, we define the differential 


operator 


é O 0 0 
f= ii ga gag Ba) 


in the natural way. In [12], Cooper and Dutle use the aforementioned Morozov- 


Shakirov formula to show that the d-th trace of Ay, 


Tra(An) =(k-D" agape] (3.1) 


where tr(A%*-)) is the standard matrix trace (for a more detailed explanation, see 
[12]). We prove our main theorem with the aid of the following reformulation of 
Equation 3.1: 

Tra(An) = (k-1)" D0) Cu(#H CH). 


HeHg 


3.2. THE ASSOCIATED DIGRAPH OF AN OPERATOR 


Recall the d-th trace of Ay from Equation 3.1, 
n hi 
Pagosa os a (i (ik = Dy wate) 
where tr Aer) is the standard matrix operation. Let Fainds ge be an addend of 
LE fe in Tra(Az). When the context is clear we suppress the subscript and simply 


write a Given @ = (41, t2,..+, ta(e—1)) let 


Ag ‘= Ai, jin Ain ig a Aiac 1)—1:td(k ake 1),i1 


and recall that 
tr(AS%P-Y) = 5 Ag 


a 


where the factors of A, are commutative. Adhering to the terminology of [12] we say 
Aq is k-valent if k divides the number of times 72 occurs in a subscript of Ag. We 
utilize divisibility notation for monomials in tr(A“*~), e.g., using g|h to denote that 


g occurs as a factor of the formal product h. We say that Ag survives f if As #0). 


Definition 3.4. For a differential operator fiechs the associated digraph of f , de- 
noted D, is the directed multigraph where there are d; distinguishable edges directed 


f and isolated vertices are ignored. 


dj 
from 2 to 7 given (a2) 
UJ 


We suppress the subscript and write D when 7 is understood. We recall the 


following graph theoretic definitions according to [14]. 


Definition 3.5. A walk in a graph G is a non-empty alternating sequence 
UVpEQU1e, ..- Ex—1Uk 


of vertices and edges in G' such that e; = {v;,v;41} for all i < k. A walk is closed if 
Up = Up. A closed walk in a graph is an Euler tour if it traverses every edge of the 
graph exactly once. An Euler circuit is an Euler tour up to cyclic permutation of its 


edges, i.e., an Euler tour with no distinguished beginning. 


We denote the set of Euler tours of a graph G which begin at the edge e € E(G) 
by €.(G) and we denote the set of Euler circuits of G by €(G). Recall that a digraph 
D has an Euler circuit if and only if deg*(v) = deg (v) for all v € V(D) and D is 


weakly connected. 


Definition 3.6. Let w = (v;)79 be a sequence of (not necessarily distinct) vertices 
of D. We say that w describes an Euler tour in D if there exist distinct edges 
€0,-++,€m—1 Such that vpeovyey...€m—1Um is an Euler tour in D. Moreover we say 


that such Euler tours are described by w. 


Note that the use of Euler tour in the previous definition is well-founded as eo is 


distinguished as the first edge. 


Lemma 3.7. Consider Trg(H), ftr(A¢%*-)) 4 0 if and only if D is Eulerian. In 
this case 


ftr(A%) = |B(D)||€(D)]. (3.3) 


Proof. Consider Trg(H). Fix a term A, of tr(A@-)) and a differential operator f 
of [Tj=4 fies Suppose fAa # 0. Whence fas # 0 the factors of 3 are in one-to-one 
correspondence with the factors of A,. It follows that the edges of D fare in one-to-one 
correspondence with the factors of A,. Notice that for i € V(D), deg* (i) = deg (i) 


by Equation 3.2. Further, D is strongly connected as the sequence of indices from 


10 


i to 7 (cyclically, if necessary) is a walk from vertex i to vertex j7. Therefore, D is 
Eulerian. 

Suppose now that D is Eulerian and let a = (v;);<9 describe an Euler tour in D. 
We claim that fA is equal to the number of Euler tours in D f described by a. Let 


m(i,j) denote the number of edges from i to j in D. Notice that 


Ab a I] Ares) and f = II i 
ate ig OA 
i,j€V (D) i,jEV(D) a,j 
Clearly 
fAc= JE mi). 
i,jeV(D) 


Moreover, a describes |]; ;ev(p) m(i, 7)! Euler tours in D by straightforward enumer- 
ation. Observe that there are |F(D)||€(D)| Euler tours in D as we may distinguish 


any edge from D as the first edge of an Euler circuit. Thus, 


Ftr(Ate) = Li Ae (D)||e(P)| 


as every Euler tour is described by exactly one a. 


We conclude this section with a remark about the evaluation of Equation 3.3. 
Conveniently, |€(D)| can be computed using the BEST theorem, originally appearing 


in [1] as a variation of a result of [38]. 


Theorem 3.8. (BEST Theorem) The number of Euler circuits in a connected Eule- 
rian graph G is 


|€(G)| = 7(G) T] (deg(v) — 1)! 


vEV 


where T(G) is the number of arborescences (i.e., the number of rooted subtrees of G 


with a specified root). 


For simplicity we abbreviate r(f) = 7(D;). Combining the BEST theorem and 


the observation that |E(D)| = d(k — 1) yields the following. 


11 


Corollary 3.9. 
ftr(A%*) =d(k—1)r(f) J] (deg-(v) - 1)! 
veEV(D) 
As a final note, recall that 7(G) can be computed using the Matrix Tree Theorem, 
which makes the computation of the right-hand side of the equality in Corollary 3.9 


efficient. 


Theorem 3.10. (Matrix Tree Theorem/Kirchhoff’s Theorem) For a given connected 
graph G with n labeled vertices, let Ay, A2,..-An—-1 be the nonzero eigenvalues of 


L(G) = D(G) — A(G). Then 


3.2.1 EULER OPERATORS AND VEBLEN HYPERGRAPHS 
In Lemma 3.7 we showed that the only differential operators f for which 
ftr(A“) 40 


are the operators whose associated digraphs are Eulerian. The question remains: 
which f have an Eulerian associated digraph? We answer this question with the 


following graph decoration. 
Definition 3.11. We define the u-rooted directed star of a k-uniform edge e to be 
Se(u) = (e,{uv:v €e,u F v}). 


A rooting of a k-graph H is an ordering R = (S_,(v1), Se,(V2),---; Sem(Um)) such 
that E(H) = {e1,...,@m} and vu; < vi41. Given a rooting of H we define the rooted 
multi-digraph of R to be 

Dp= U Se,(vi) 


i=1 


where the union sums edge multiplicities. We say that a rooting R is an Euler rooting 


if Dr is Eulerian. We denote the mutli-set of rooted digraphs of H as S(H). 


12 


Note that two distinct rootings can yield the same rooted digraph. We suppress 
the subscript Dr and write D when the context is clear. We further refer to D as a 


rooted digraph of H for convenience. 


Definition 3.12. Given a rooted digraph D € S(H), we define the rooted operator 


of D to be 
* O 
La OAuy 


uv€ E(D) 


Moreover, we denote 


S(H) = {fo : D € R(H)}. 


In the case when D is Eulerian we refer to fee as an Euler operator. 


The notation of fp is consistent with our usage of - whence 
fo\ IF 
i=1 
where d; is equal to the number of times vertex i is appears as a root of D. If f is 
a rooted operator then it is understood that there exists a (not necessarily unique) 


rooting R such that, with a slight abuse of notation, f = fis We call such a rooting 


an underlying rooting of a differential operator. 


Lemma 3.13. The associated digraph of an operator fg is Eulerian if and only if f 


is an Euler operator. 


Remark 3.14. By Lemma 3.13 the only operators which have nonzero contribution 
to Trg(H) are Euler operators. We denote R(H) C S(H) to be the multi-set of Euler 


rooted digraphs of H. We further denote 
R(H) = {fo : D € R(H)}. 


Remark 3.15. One can deduce Theorem 4.1 of [35] from Lemma 3.13 by a change 
of notation: our fp is their F', our set of Eulerian associated digraphs arising from 


fris their Egm_1(n), and our €(D) is their W(E). 


13 


We now show that an Euler rooting is a rooting of a special type of hypergraph. 
Definition 3.16. A Veblen hypergraph’ is a k-uniform, k-valent multi-hypergraph. 


Lemma 3.17. An Euler rooting R is a rooting of precisely one labeled Veblen hyper- 


graph. 
Proof. Suppose S = (S;)'@, is a rooting of a connected k-graph H. Since Dg is 
Eulerian we have for all 7 € V(H) 
deg*(j) = (k— Di fi: v1 = Fh] = is ws A 5,9 € ex}| = deg (4). 
Fix a vertex v € V(H). We compute 
degy(0) = degh(v) + deep (v) 
=|{i:v,=v}|4+ |{t: 0, 4 v,v € e;}| 


=|{i:v,=v}|+ (k-LD|{i: v; = v}| 


= hifi ts= vp). 


Observe that k | deg,,(v); it follows that H is Veblen by definition. Now suppose 
that Ho is a connected Veblen graph such that S is an Euler rooting of Hp. As 


S is a rooting of Ho, E(Ho) = E(H) and since both hypergraphs are connected 


V(Ho) = V(H). It follows that H is unique. 


We combine Lemmas 3.13 and 3.17 into the following Lemma. 


Lemma 3.18. We have 
fir(A"™) £0 


if and only if f = fie is a rooted operator. Moreover, the underlying rooting of f is 


necessarily an Euler rooting of precisely one connected, labeled Veblen hypergraph. 


‘The nomenclature is a reference to Oswald Veblen (1880-1960) who proved an extension of 
Euler’s theorem in 1912. We present a brief note about Veblen’s namesake theorem at the conclusion 
of this section. 
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In the following section we use Lemma 3.18 to express the codegree-d coefficient of 
a k-graph as a function of Veblen hypergraphs. Here we conclude with a note about 


Veblen’s theorem. 


Theorem 3.19. (Veblen’s theorem [39]) The set of edges of a finite graph can be 


written as a union of disjoint simple cycles if and only if every vertex has even degree. 


Unfortunately, Veblen’s theorem does not extend to higher uniformity: the set of 
edges of a finite k-graph H can not always be written as a union of disjoint simple 


k-regular k-graphs if and only if H is k-valent. Consider the Veblen 3-graph 7 which 


consists of three bottomless tetrahedrons each sharing a common base. To be precise, 


T= ((o b,¢, 1.2.3} U (ea waar? 3) . 


A drawing of 7 is given in Figure 3.1. Since there are only three edges containing 


2 


Figure 3.1 The Veblen 3-graph 7 where edges are drawn as triangular faces. 


i € [3] it must be the case that any partition into Veblen graphs places each edge 
containing 2 into the same class. Observe that for each 7 the vertices a, b, and c each 


have degree 2. Therefore, the only 3-valent edge partition is the trivial one. 


3.3 THE ASSOCIATED COEFFICIENT OF A VEBLEN HYPERGRAPH 


We now turn our attention to computing the codegree-d coefficient of a k-graph H 


via Equation 3.1. From Lemma 3.18 we know that the only operators which satisfy 


15 


£ tr(A“*-)) 4 0 are rooted operators. Furthermore, as a differential operator of 
Tra(Axz) is of degree d, the underlying Euler rooting of ic is a rooting of precisely one 
connected, labeled Veblen hypergraph with d edges. We equate Trg(H) to a weighted 
sum over Euler rootings of connected Veblen graphs with d edges which “appear” in 


H. Consider the following generalization of the notion of subgraph. 


Definition 3.20. For a labeled multi-hypergraph H, we call the simple k-graph 
formed by removing duplicate edges of H the flattening of H and denote it H. We 
say that H is an infragraph of H if H C H. Let Va(H) denote the set of isomorphism 


classes of connected, labeled Veblen infragraphs with d edges of H. 


Definition 3.21. The associated coefficient of a connected Veblen hypergraph H is 


TD 
Cy = : 
‘ ok los. ) deg™ (v) 


The associated coefficient of a (possibly disconnected) Veblen hypergraph is 


CH — [| Ca. for H = LJ Gi. 


i=1 i=1 
We present the associated coefficient of Veblen 3-graphs with six or fewer edges 


in Table 3.3. 
Definition 3.22. For a k-graph H and a Veblen k-graph H = Uj, Gi; we define 


(#H CH) = LL lana Olt e Ha =| Gall 


‘Ant(E)| 


In the case when H is connected this simplifies to 


Ger cy — luli) 


SH) = Rants SM: 8 = BH = [Aut(H)/Aut(#)|- {8 SH: = HB, 


Note that for H = U2, G;, (H C H) is not multiplicative over the components of 


H as 


us _ (#H CH)|Aut(H)| 
Il4G: ¢ ) = TTT OTTGR Ta 


i=1 


However, we have the following identity. 
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Table 3.1 The associated coefficient of Veblen 3-graphs with six or fewer edges. 


B(H) GH 
[0, 1, 2], [0, 1, 2], (0, 1, 2] 3/8 
[0, 1, 2], (0, 1, 3], (0, 2, 3], (1, 2, 3] 21/8 
[0, 1, 2], [0, 1, 4], (0, 3, 4], (1, 2, 3], [2, 3, 4] 51/16 
[0, 1, 2], (0, 3, 4], (0, 3, 4], (1, 2, 3], [1, 2, 4] 27/16 


0, 1, 2], (0, 1, 2], [0, 1, 2], (0, 1, 2], [0, 1,2], (0, 1,2) | 3/16 

0, 1,21, (0° 1, 2), [0; 1,2), 10, 1,3), (0,1, 3],10;1;3) | .9/8 

0, 1, 2], (0, 1, 2], (0, 1, 2], (0, 3, 4], [0, 3,4], (0, 3,4) | 9/32 

0, 1, 2], (0, 1, 2], (0, 1, 3], [0, 2, 4], [0, 3, 4], (0, 3,4] | 99/32 
0, 1, 2], (0, 1, 3], (0, 1, 4], [0, 2, 3], (0, 2, 4], (0, 3, 4] | 213/16 
0, 1, 2], (0, 1,3], (0, 4, 5], [1, 4, 5], [2, 3, 4], [2,3,5) | 69/16 
0, 1, 2], (0, 1, 3], (0, 3, 4], [1, 3, 5), [2, 4, 5], (2,4,5)] | 63/32 
0, 1, 2], (0, 2, 3), (0, 3, 4], [1, 3, 5], (1, 4,5], [2, 4,5] | 129/32 
0, 1, 2], (0, 1, 2], (0, 1, 3), [2, 4, 5), [3, 4, 5], (3, 4,5) | 27/32 
0, 1, 2], (0, 1, 3], (0, 2, 3], [1, 4, 5], [2, 4, 5], (3, 4,5] | 63/16 
0, 1,2], (0, 1, 3], (0, 2, 4], [1, 3, 5], [2, 4,5], [3, 4,5] | 117/32 


Lemma 3.23. Let H = U2, Gi be a Veblen k-graph. If uy denotes the number of 
linear orderings of the components of H (where two components are indistinguishable 


if they are isomorphic) then 


(HH CH) ="* THE CH). 


Os 
Proof. Suppose there are t isomorphism classes of components of H with represen- 
tatives H,, Ho,...H;,. Denote the number of components of H which are isomor- 
phic to H; as p;. Fix an ordering of the components which are isomorphic to H;, 
{Gi i,Goi,-.-,Gy,i}. The number of distinct linear orderings of the components of 


H where G‘; and G; are indistinguishable when G,; = G,,; is 


( m 
MH = 
L41, L2, cay Lt 


m)\ : 
— = al. 
ya LP 


so that 


Le 


Note that for a € Aut(#), there exists o € G,, such that a(Gj,) = Go.;);- In this 
way, Aut(H) induces a permutation on the isomorphism classes of the components 
of H (note that this map is well-defined since the components are labeled). Let 
w: Aut(H) > G,, x Gz. X ++: x G,, be such a map. Notice ker(w) is the group of 


automorphisms of H which maps each component to itself. Appealing to the First 


Isomorphism Theorem we have 


Aut(H)) ate 
Ti, aut(@y) ~ LL 


The desired equality follows by substitution. 


Remark 3.24. The equation in Lemma 3.23 implies that (#H C H) is multiplicative 


over its components if and only if the components of H are pairwise non-isomorphic. 


Let 
A(d,n) = {(di,..-,dn) : od: = d,d; > 0} 


be the set of arrangements of d into n non-negative parts and further let A*(d,n) C 
A(d,n) be the set of arrangements of d into n positive parts. For a k-graph H and 
a € At(d,|E|) let R*(H) be the set of Euler rootings of all labeled, connected Veblen 
infragraphs of H which have the property that vertex vu; is the root of exactly d; 
edges. (N.B. We take a € A*(d,|E]) as it is necessary that d; > 0 for D € R*(H) to 


be Eulerian.) 


Remark 3.25. Let Vi(H) denote the set of (possibly disconnected) Veblen infra- 
graphs of H with d edges up to isomorphism. Further let Va(H) C Vi(H) denote the 


set of connected Veblen infragraphs of H with d edges up to isomorphism. 


We now present a formula for Trg(H) as a weighted sum over its Veblen infra- 


graphs. 
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Lemma 3.26. For a k-graph H 


Tra(H) =d(k—-1)" 3 Cu(#H CH). 


HEVa(H) 


Proof. For convenience let |E'| = |E(H)|. We equate 


TD 
S cu#Hcn= > ( 3 : ) enc) 
Heva(H) Hev4(H) \ \DeR(H) Twevpy) deg” (v) 
TD 
- ¥ s : ) 
a€A+(d,|E|) \DER*(H) Teva) deg” (v) 
Recall 


Tra) =(k-Y (Hage), 


Applying Lemma 3.18 we have 


Tra(H) =(kK-1)"" SP 


ac A* (d,|E}) 


fp inl Ae 1) 
pereny Wei (4i(® — 1)! 


By Corollary 3.9 , 


fotr(A%)) =d(k-1)tp J] (deg @(v) — 1)! 
ve V(DrR) 


When D € R°(H) with a = (d,,...,dn), we have deg” (v;) = d;(k — 1). By substitu- 


tion we have 


Tra(H) = d(k — 1)" 


= (5 maces) 
acAt(ajE) \wer@(H) [lvev(ny des” (v) 


=d(k—1)" S$) Cu(#H CH). 
HEVa(H) 


We are now prove a generalization of the Harary-Sachs formula for k-graphs. 


Theorem 3.27. For a simple k-graph H, 


d(H) = Y> (-(K- 1)")Cu(#H CH). 


Hevi(H) 
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Proof. Fix a k-graph H and d > 1. From [12] we have by Equation 3.1 


bal) = Py (= Tri(H) — Tro(H) — ae) 


i : 2 re d 


where 
d 
Pi(ti,t,--..t => YS aa 


m=l dit -+dm=d ae 
Fix an 1 < m < d and an arrangement a = (d,,do,...,dm) € At(d,m). By 
Lemma 3.26, 


= Tra, (H) 


oe =—(kK-1)" YY Col#G CH). 


GEeVa; (H) 
Let V*(H) denote the set of m-tuples of connected unlabeled Veblen infragraphs of 
H whose i-th coordinate has d; edges for i € |m], such that 30; d; = d. We have 


fe = (yey TT cower cH) 
i=1 4 i=1 Geva,; (H) 
=(-(k-I" Cn TG CH) 
H=G1U:-UGm i=l 


(G1,G2 seess Gm)€V*(H) 
For m € N, let V"(H) be the set of unlabeled Veblen infragraphs of H with d 


edges and m components. Appealing to Lemma 3.23 we may write 


a Tr(H)  Tre(H) Trg(H) 
a= Fa(- 1 ae 2 vo 
ak 1 -Tra(H) 
oe enw LLG 
= > (-(k-1))" Yo Cy (Tei (#Gi “)) 
ace At(m,d) HEV2(H) m\ 


Me iMs 


3 
I 


Coe cw (Ht TTeec < %)) 
Peace 


Hey” 


= (—(k — 1)")Cu(#H CH). 
Hev*(H) 
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3.4 DEDUCING THE HARARY-SACHS THEOREM FOR GRAPHS 


With Theorem 3.27 in hand we can express the codegree-d coefficient of the normal- 
ized adjacency characteristic polynomial of a hypergraph as a weighted sum of Veblen 
infragraphs with d edges. In the case when G is a 2-graph our theorem simplifies to 
galG) = So (-1) Cu (#H CG). 
Hev3 
Recall that the Harary-Sachs theorem (Theorem 3.1) expresses the codegree-d coeffi- 
cient as a weighted sum over certain subgraphs on d vertices, whereas Theorem 3.27 
expresses the same quantity as a weighted sum over certain subgraphs with d edges. 
We now argue that these two sums are equal. 

An elementary subgraph of a graph G, is a simple subgraph of G whose compo- 
nents are edges or cycles (see [3] for further details). In keeping with their notation, 
let Ag(G) be the set of elementary subgraphs of G with d vertices. Notice that a con- 
nected elementary graph is the flattening of a cycle (e.g., the flattening of a 2-cycle is 
an edge). Recall that cycles (and disjoint unions of cycles) are the only two regular 


non-empty graphs which have an equal number of vertices and edges. Indeed 
Ag(G) C {Hs HE Vi(@)}. 


By straightforward computation we have that the associated coefficient of a 2-cycle 
is 1 and the associated coefficient of a simple cycle (i.e., any cycle which is not a 
2-cycle) is 2. Restricting our attention to Ag(G), we have by Theorem 3.27 
(1) Ca#H € G) = (—1) 2) (HH CG) (3.4) 
HeEAg(G) 
where z(#) is the number of cycles in H. Note that Equation 3.4 is the conclusion 
of the Harary-Sachs theorem. We deduce the Harary-Sachs theorem from Theorem 


3.27 by showing that the summands of 


pal 


which do not arise from elementary graphs sum to zero. We make this statement 


precise with the following. 


Definition 3.28. For a multigraph G and an edge e € E(G), write m(e) = me¢(e) 
for the multiplicity of e in E(G). Let G be a connected, labeled Veblen graph with 
distinguishable multi-edges. Given a multiset P of multigraphs whose multi-edges are 
indistinguishable, each on the vertex set V(G), we write PF G if, for each e € E(G) 
S> mp(e) = me(e). 
P,EP 


Lemma 3.29. 
1 :G is a 2-cycle 
3 (—1) Cp = 4 2 :G is a simple cycle 
P-HG 
0 : otherwise. 
Remark 3.30. We refer to Veblen 2-graphs as Veblen graphs. An Fuler orientation 
of a graph G is an orientation of the edges of G such that the resulting digraph is 


Eulerian. 


We first provide a combinatorial formula for the associated coefficient of a Veblen 


graph. 


Lemma 3.31. Let G be a connected Veblen graph. We have 


|E(@)| 


Cg = =. 
Ileex(a) mie)! 
Proof. Let G be a Veblen graph. Since G is Eulerian, we write deg(v) = 2d, for 


convenience. Because deg (v) = d, for all D € R(G) we have 


i) 
Ce SSS S- TDs 


Ivev(a) dy DER(G) 
By the BEST theorem (i.e, Theorem 3.8) we have 


_ ey 
Teva (de a 1)! 


TD 


2 


Let Np(v) denote the out-neighborhood of v in D and let degp(v,u) denote the 


number of edges directed from v to u in D. We denote 


dy a! 
Gy 2 Tluenp(v) degp(v, u)! 
which is the number of linear orderings of out-edges of v in an Euler orientation D. 
Consider the equivalence relation ~gg) where R ~ R’ if and only if De = Dp. 
Note that ~ identifies two Euler rootings if their associated digraphs are the same 
Euler orientation. Let [R] denote the equivalence class of R under ~. Suppose 


R(G) = Us_,[Ri] and note 


le ie) 


as two rootings in [A] differ only in the ordering of the u-rooted stars for u € V(G). 
Let O(G) = {D,,..., Di}, where D; = Dpr,, denote the Euler orientations of G. For 


convenience, we write N; for Np,. We equate 


t dy 
vere _ Xl a “ (: veV(G) Ga) | 


Substitution and simplification yields 


Je(D,) | 
i=1 [vevia) (Tuer) deg; (v, u)!) 


Ce = 


Since deg,(u, v) + deg;(v, u) = m(uv) we have 


( (wv) )- (wv)! 
deg;(u,v))/ deg; (u, v)! deg; (v, u)! 
so that 


_ e)! 


m( 
Iheviay Tuevia) deg;(v, u)! Tucvev(G eG u, v)! deg;(v, u)! 
v) 


uve€ E(G),u<v Con (u, v) ) 


Ilee x(a) m(e)! Tleex(a 
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Then 


t 
|E(D)| 
C= 
‘ » Iloev(a) (TMuens( deg; (v, w)!) 
7 Teena) me)! 3 |E(Di)| 


Teema) Me)! 4 Teva) (Tucncts deg;(v, u)!) 


hese Wye * Mmg(uv) | 
7 Heee@) ane) 2 (_. — ee 
|E(@)| 


Ileex(a) mie)! 


where the last equality follows from the observation that D; has indistinguishable 


multi-edges. 


We now prove Lemma 3.29. 


Proof. Let G be a connected Veblen graph which is not a cycle. Further assume that 


the multi-edges of G are distinguishable. We aim to show 


S~ (-1)°Cp = 0. 


PFG 


By Lemma 3.31 we have for connected G, 


|E(G)| 
Ce = =". 
= Tleex(ay m(e)! 


Let P = Ui_, P; be a disjoint union of Veblen graphs. We denote 


ea & me(e)! 
P(e) J Tha1mep,(e)! 
and 

le(P)| = Tle). 


We equate 


II me S>(-1)")Cp = Yo (-1) €(P)|_ T] ee 


e€E(G) 
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Notice that |€(P)| T]eexa) () counts the number of partitions of E(G) into edge- 
disjoint Euler circuits of graphs on V(G) with unlabelled edges which are precisely 
the elements of P. 

We say that an Euler circuit is decomposable if it can be written as a union of (at 
least) two edge disjoint Euler circuits, and indecomposable otherwise. The number of 
decompositions of Euler circuits of G into exactly t parts is 

» (ew I oy | 
PEG, e€ E(G) 


By Inclusion/Exclusion, the number of indecomposable Euler circuits of G is 


iy (ie I ) , 
t=1 PEG, e€E(G) 

We assumed that G is a Veblen graph which is not a cycle. We have by Veblen’s 
theorem (i.e., Theorem 3.19) that every Euler circuit in G is decomposable. It follows 
that G has no indecomposable Euler circuits; that is to say, 


SP Ds (ie II =) =0, 


t=1 PEG, e€E(G) 


from which the desired conclusion follows. 
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CHAPTER 4 


LOW-CODEGREE COEFFICIENTS 


From the Harary-Sachs Theorem (i.e., Theorem 3.1) we know that for a 2-graph G, 
if O(G) = +9 ¢2"' then 


Co = 1,¢, = 0, co = —|E(G)|, and c3 = —2(# of triangles). 


The goal of this chapter is to provide a similar description for the low-codegree coef- 
ficients of a k-graph. We begin by providing such a description for the first & proper 
coefficients of d(H). The following is a collection of results from [12] but we provide 


a new proof via Theorem 3.27. 

Lemma 4.1. For a k-graph H and d € |k — 1], 6(H)q =0. Moreover 
Ox(H) = —(k — 1)" *k* | E(H). 

Proof. Recall from Theorem 3.27 we have 


d(H) = SY) (-(k- 1)")Cu(#H CH). 
Hevi(H) 


Notice that for d € [k—1], Vi = @ so that ¢a(H) = 0. Further, notice that Vi = {Ex} 


where EF; is the k-uniform edge. By Lemma 4.2 we have that 


R(Ex) = {((1, [A]), (2; (AI), ++ (Rs (AI) F- 


For R € R(E;) we have Dp ~ K;, so that L(Dr) = kl, — J, whose multi-set spectrum 
is {0,k,k,...,k}. By Kirchoff’s Theorem we have Tg = k*~? so that 


kk-2 
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Substitution yields 


ch = (—(k — 1)")Cu, (#ER CH) = —(k — 1) “kK E(H)]. 


We begin by providing insight into the computation of Cy. Then we turn our 
attention to an explicit formula for ¢,41(H). In the final section we provide an explicit 


combinatorial formula for the first six proper leading coefficients of a 3-graph. 


4.1 COMPUTATIONAL NOTES 


For an arbitrary Veblen k-graph we compute C'y by first determining a list of all 
Euler rootings. Then for each Euler rooting we compute the corresponding summand 
via Kirchoff’s Theorem. The difficulty in this computation is in determining the set 
of all Euler rootings efficiently. The following lemma imposes a necessary condition 


on rootings. 


Lemma 4.2. Let H be a Veblen k-graph. Fixt € N andd such that tk < d < (t+1)k. 


Let fie a a, be a differential operator of Tra(H). If d; >t then 
f tr(A%@-Y) =0. 
In particular, if f tr(At@-D) #0 then d; < Fal and this bound is sharp. 


Proof. Consider Tra(H) and fix a differential operator fue d,. Observe that for 


D = D; we have deg} (i) = (k—1)d; and deg, (i) < d—d;. Without loss of generality, 


suppose d, >t. Then 
deg (1) <d—d, <d—(t4+1) < (t+ 1)k-(¢4+1) = (¢4+ 1)(k—-1) < deg$(1). 
By Euler’s theorem D does not have an Euler circuit. Appealing to Lemma 3.7, 
f tr(A%*-)) =0 


yas 


and the first statement follows. 

we have by our first statement d; < t < |d/k|. To see that this bound is sharp, fix 
t€ Nand k > 2. Set d:= kt and consider Trg(H) where H = (kj, {d x [k]}) is the 
k-uniform edge with multiplicity d. Choosing d; = t for each i € [k] yields precisely 


one differential operator f as 


| 


vj 
Observe that D is the complete multi-digraph on k vertices where each pair of vertices 


has t edges oriented in both directions. By Euler’s theorem D has an Euler circuit. 


Let a@ describe such a circuit. By our previous claim A, survives f Sie) f tr( AMD) 40 


as desired. 


Note that, for the classical Harary-Sachs Theorem, the number of elementary 
graphs (explained in the next section) one needs to sum over for the codegree-d 
coefficient is simply the number of partitions of d into positive parts. The number of 


Veblen 3-graphs is exponentially larger. 


Remark 4.3. Lemma 4.2 is most useful when d is close to k in value. In particular, 


if d < 2k then d; = 1 for all 7. 


Recall that Vz and V3; denote the number of connected and (possibly) disconnected 


Veblen 3-graphs with d edges. We have computed 
(MINS = 0,0, 172,11, 296: 1997815...) 


and further 


(\V5|)22, = (0,0, 1, 1, 2, 12, 27, 125, 795,...) 


see A320648 [36]. 
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4.2. THE CODEGREE-(k + 1) COEFFICIENT OF A k-GRAPH 


nti (en(tS") 


be the k-uniform simplex. It was shown in [12] that the k-uniform simplex is the only 


Fix k > 2 and let 


connected Veblen hypergraph with k +1 edges, up to isomorphism. Moreover, it was 


shown for a k-graph H 
bin (KO,) = (hk - 1)" OH CH) 


for some constant Ci; depending only on k. The authors of [12] were able to show 
that Cp = 2, C3 = 21, Cy = 588, C's = 28230 via laborious use of resultants. In this 
section we provide an explicit, efficient formula for C, and use it to compute Co. 
To that end, we first describe the set of rootings of the k-uniform simplex. 


For the remainder of this section let H denote the k-uniform simplex Iv ae 


x @k-1 


Observe that the 7-th coordinate polynomial of Ay can be written 


UX... UK41 
a a 
ifi i 5 
fori € [k + 1]. Let Dy41 be the set of derangements of [k + 1], i-e., permutations 


without any fixed points. Recall that S.(u) is the u-rooted directed star of e. 


Lemma 4.4. Foro € Dz. define 


k+1 
De = VU Stes gowm} 
i=l 


Then 
R(A) = {Do :0 € Dri}. 


Proof. Let o € Dg11 and consider D,. We suppress the subscript and write D for 


convenience. Note that D is Eulerian whence 
degp(j) = (k—1)|{i: 3 = j}| =k -—1=|{i: 4 45,9 € es}| = degp(s) 
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for all 7 € [k+1] as there are exactly & edges which contain j and o is a derangement. 
Suppose D € R(H). As D is Eulerian we have by Lemma 4.2 
k+1 
D= VU Sterno. 


i=1 


In particular, {v1,...,ve41} = [k +1]. Observe that o = {(i,u,)}E) © Derr is a 


derangement. The conclusion follows from the fact that D = D,. 


This immediately implies the following. 


Lemma 4.5. For the k-uniform simplex H, 


To 


Cy= ——— 
o€Dx,, Llvev(p.) deg (v) 


We now show that a summand in the aformentioned formula of Cy depends only 


on the cycle type of the derangement. 


Theorem 4.6. Let o € Dzyi1 with cycle decomposition c,c2...c, where cycle c; has 


length £;. Then 


1 T (pei e+1 
Ca= GI DFI(R+D) S- II (« + (-1) Ms 


O=C{CQ...Ct Dr+1 i=l 


We first prove a technical lemma. The notation spec(/) denotes the ordinary 


(multiset) spectrum of a matrix M. 


Lemma 4.7. Foro € Sp41, spec(M, — J) = (spec(M,) \ {1}) U{—n} where M, is 


the permutation matrix associated with o. 


Proof. Let o be a permutation of [n+1] with cycles c,, co,...,c; of length ly, lo,...,, 


respectively. Recall that the spectrum of M, is given by 


spec(M,) = Lic, Gr ae et 


i=1 


30 


and note that the spectrum of M, depends only on the cycle type of 0. Without loss of 
generality, suppose that the cycles of o are increasing (i-e., [1, 1], [41 +1, 02+4],...). 


Consider the following block partition 


By, —J, —Jo, 

—Jp, Bo —Je, 
M,- n+1 = 

95, Se, ihe. Be 


where B; is the 1; x 1; square circulant matrix corresponding to «;, 


-1 0 -1 -1 
-1 -1 0 -1 
B= 
-1 -1 -1... 0 
Q -1 -1... -1 


Let 7 > 0 and consider the eigenpair (G x) of M,. We compute 
(M, — J)x = M,x —- Jx = Gx —-0= Gx 


where Jx = 0 because the coordinates of x corresponding to C; are the complete set 
of @;-th roots of unity. 
We now show 1 € spec(M, — J) has a geometric multiplicity of t—1. Fix 1<i< 


t — 1 and consider x € C”*! where 


1 :7 EC; 
oe ai =a ye Cag 
0 : otherwise. 


For j € C; we have, 


é, 
7 (Gti) = 1= xy, 


(Me ~ Dx) = “(GY + 


dl 


for 7 € C;41 we have, 


b; e; 
(Me — J)x)j = 4 + 7 (Gin — YI = FP = 
i+] i+] 
and for 7 ¢ Ci, C41 we have 
f, 
((M. — J)x); = -4 + 7—(Gin1) = 0 = xj. 
fa 


Therefore, (1,x) is an eigenpair for 1 < i < t; moreover, these vectors are linearly 
independent. 


Finally, consider the all-ones vector 1 € C"*! where 
(M, —J)1l1=-—nl 
so (—n,1) is an eigenpair of M@, — J. We have shown 
spec(M, — J) 2 (spec(M,) \ {1}) U{—n} 


and the reverse inclusion follows from the observation that the multiplicities on the 


right-hand side add up to (at least, and therefore exactly) n + 1. 


We now prove Theorem 4.6. 


Proof. Consider 


ges ks 


o€D141 Lvev(De) deg (v) 


Applying Lemma 4.2 (choosing d = k + 1) implies that each vertex of H is a root of 
exactly one edge in an Euler rooting of H. It follows that, for all o © Dp41, 


I] deg (v) = (k-1)8"". 


veEV(Dao) 


Thus, 
Cpa Se. 


= k+1 
Cs 1) - TCD R41 


Consider o = c,c2...c; where cycle c; has length @;. Observe that 
L(Ds) =kI+M, — J. 
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Since spec(kI + M, — J) =k+spec(M, — J) we have by Lemma 4.7 and Kirchoff’s 


theorem 
Mas (HY + (-1)*") 
k+1 


i= 


The desired equality follows by substitution. 


Remark 4.8. It was shown in [12] that 

bryi(H) = —C,(k — 1)"* 
and by Theorem 3.27 we have 

bryi(H) = —Cy(k — 1)”. 
In our notation, we can write Cy = (k — 1)*Cy. 


For ease of computation we consider Ci instead of Cy. As it is stated, Theorem 
4.6 is slow to compute as |M,41| ~ n!/e. However, summing over all derangements 
of [k + 1] is wasteful, as we have shown C;, is a function only of the cycle structure of 


a. In fact, 
t 


[eee ) 


i=1 


is constant for derangements with the same cycle type. We present a reformulation 


of Theorem 4.6 which has the advantage of considering a smaller search space. 


Definition 4.9. Let P(n) be the set of partitions of n and let Ps2(n) C P(n) be 
the set of partitions of n into parts of size at least 2. For p € Pso(k +1), let 
Drii(p) C Dei be the set of derangements whose cycle lengths agree with the 
parts of p. Further, for a partition p € P(n) let V, : [n] > [0,n] be the map 


Vet) = 19 + Dy = OF. 


We reformulate Theorem 4.6 (for C;) as follows. 
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Corollary 4.10. 


» (is) ll (bP + (-1)"*1)) 


) p=(p1,-.pt)€Pso(k+1) i=1 


where 
(k +1)! 
Tis pi) (TIE Vp (6) 


[Di+1(p)| = 


Proof. It is sufficient to show 


(k +1)! 
Ti, ps) (TIE? Vp (3)! 


for p = (pi,.--, pt) € Poo(k +1). Let A: Gy41 > Deyil(p) by 


[De+1(p)| = 


A(o) = (o(1),0(2),--.,0(p1))(o(pr + 1), ---,0(p1 + Bo) 


Note that A is surjective. Since a cycle of A(ca) can be written with any one of its p; 
elements first and there are [[**} V,(z)! linear orderings of the cycles by non-increasing 


length we have for 6 € Dri4 


t k+l 
JA“) = (It) (I vita) 
i=1 i=2 
In particular, |A~'(6)| is constant for 6 € D,41(p) so that 


ISenl= D7 [A“*()| 


6€Dx41(p) 
which implies 
Siva] 
» == : 
| +1 (P)| |A-#(6)| 


Remark 4.11. Corollary 4.10 reduces the number of summands in the computation 


of C, exponentially because 
log |Ps(a)(n)| < log |P(n)| =m y/2n/3, 
but log |Dn41| & nlogn. 
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We have computed the first few values of Cy to be 


C6 = 2092206 

Cz = 220611384 

Cg = 31373370936 
Cy = 5785037767440 


Cio = 1342136211324090 


C09 = 3433452419824795908447 767 1758634630345 268960989035871 1139013913758 
7799578881707 16788865639598053642953208929209278848309297069686374206 
618031496101898485314300253248855334075609527915686375386625810970778 
814182546067369319275314946445603388 1 15577892354872286012782651661555 
3106527369037 12206018668653541524263903668524799914172228056595466145 
2080249009900 

so that C09 © 3.433...-10°48, see A320653 [36]. Note that Cy = 2 gives the well-known 
result that, for a graph G 

$3(G) = —2(# of triangles in G). 
We conclude by presenting the asymptotics of C. 


Theorem 4.12. C, ~ (k + 1)!k**! so Cy = exp(k log k(2 + o(1)). 


Proof. As o is a derangement it follows that the length of each cycle is at least 2. 


Notice 


t 
ial < (Te +!) < Paal— 0 


o=C1C2...c¢E Dp 41 \i=l 


which implies 


» 1 ki + 1") os (kU) 


o=C1C2...c¢EDp 41 \i=l € 
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Thus, 
Cr _ Derilk*t? 1 


JUG OLerDe « 


and we have Cy ~ (k + 1)!k*T1. 


4.3. LOW-CODEGREE COEFFICIENTS OF 3-GRAPHS 


We provide an explicit formula for the first six proper leading coefficients of the 
characteristic polynomial of a 3-graph. Let H be a simple 3-graph with n vertices. 


We write (# H € H) =|{S CH: S 2H)|. From [12] we have that c; = 0, co = 0, 


cz = —3-2" 3*(#e EH), 


where e = K. 8) is the single-edge hypergraph and 


C= DD KO ea, 


We use Theorem 3.27 to situation these results and further provide an analogous 
description of cs; and cg. Clearly there are no Veblen 3-graphs with one or two edges 
as each vertex must have degree at least three. It follows that c,,co = 0. There is 
only one Veblen 3-graph with three edges, the single edge with multiplicity three. By 


Theorem 3.27, then, we have 
no 
C3 = —2 gle € H). 


Similarly, [12] observed that the only Veblen 3-graph with four edges is the complete 


graph, which implies 
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where [51 is the tight 5-cycle and [59 is the 3-pointed crown (given explicitly in 


Figure 4.1). We further compute 


_— 92n 9 (#e€H) 
eo al 2 


99 213 69 63 

rae 2(#Tes € H) + Te (#P oa € H) + se (#P os € H) + 25 (#T 66 € H) 
129 of 63 i es 

tay (#167 € H) + 35° 2(#T 68 EH) + 16 (Tos EH) + 35 (#2 6,10 eH)) 


or (See H) + HT or €H) + HT o2 €H) 


where Ig; are provided in Figure 4.1. 


Table 4.1 Some connected Veblen 3-graphs and corresponding values. 


r E(T) Cr | |Aut(P)|/|Aut(D)| 
Ts. | (123)(125) (145) (234) (345) 51/16 1 
Ts2 | (123)(145)(145)(234)(235) 27/16 1 
ea (123)3(124)3 9/8 1 
Eas (123)3(145)3 9/32 1 
oe (123)2(124) (135) (145)? 99/32 2 
Tea | (123)(124)(125)(134)(135)(145) | 213/16 1 
Tes | (123)(124) (156) (256) (345)(346) | 69/16 i 
Tee (123) (124) (145) (246)3 63/32 1 
Te7 | (123)(134)(145)(246) (256)? | 129/32 1 
ie (123)2(124) (356) (456) 27/32 2 
Teo | (123)(124)(134)(256)(356) (456) | 63/16 1 
To10 | (123)(124) (135) (246) (356) (456) | 117/16 1 
Tas (123)5(145)3 9/32 2 
Es (123)3(145)3(246)3 9/8 1 
Pag (123)3(145)3(167)3 81/128 1 
igen (123)9(145)3 9/32 2 
ae (123)5(145)5 27/64 1 
Bs (123)9(145)3(167)3 81/128 3 
Tou (123)5(145)3(246)3 63/32 3 
ioe (123)3(145)3(167)3(246)3 459/64 1 
[se (123)3(145)3(246)3(356)3 255/16 1 


Remark 4.13. The Fano Plane is a Veblen 3-graph with seven edges. The associated 
coefficient of the Fano Plane is 87/16 = 5.4375. 
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CHAPTER 5 
STABLY COMPUTING THE MULTIPLICITY OF KNOWN 


Roots 


Given a k-graph H we aim to compute ¢(H). Unfortunately, the resultant is known 
to be NP-hard to compute [22]. An efficient method for computing the resultant, even 
in special cases, would impact several fields of mathematics, perhaps nowhere more so 
than computational algebraic geometry. Performing this computation directly using 
built-in routines has proven to be intractable even for 3-uniform hypergraphs with few 
vertices. Nonetheless, one can attempt to imitate classical approaches to computing 
characteristic polynomials of ordinary graphs. 

A notable example of this is the work of Harary [23] (and Sachs [33]) who showed 
that the coefficients of ¢(G) can be expressed as a certain weighted sum of the counts 
of subgraphs of G. This formula allows one to compute many low codegree coefficients 
—i.e., the coefficients of «4* for k small and d = deg(¢,) — by a certain linear 
combination of subgraph counts in H. Alas, this computation becomes exponentially 
harder as the codegree increases, making computation of the entire (often extremely 
high degree) characteristic polynomial impossible for all but the simplest cases. 

We could instead view the leading coefficients as the result of mixing minimal 
polynomials of ¢(H) according to their multiplicity. We make this problem precise 


as follows. 


Problem 5.1. Let K be a field of characteristic zero. Is it true that a monic poly- 


nomial p € K[zx] of degree n with exactly & distinct, known roots is determined by 
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its k proper leading coefficients? 


We answer Problem 5.1 via an algorithm which allows us to stably compute the 
multiplicity of each root exactly. The algorithm is stable in the sense that if an eigen- 
value is approximated by an ¢-disk, where ¢ depends “reasonably” on the parameters 
of the problem, the resulting disk approximating its multiplicity contains exactly one 
integer. Furthermore, the computational complexity of the algorithm is quasi-linear 
in nb where b = O(Ine) and n is the number of roots. We address the problem of 


computing o(H) in the following chapter. 


5.1 DETERMINING THE MULTIPLICITIES OF KNOWN ROOTS 


Given the set of roots of a polynomial without multiplicity and an appropriate num- 
ber of leading coefficients, one can determine the multiplicity of its roots using the 
Faddeev-LeVerrier algorithm, a matrix form of the Newton Identities. We analogize 


this result as follows. 


Theorem 5.2. Let p= 0") Ga"! € Kz] be a monic polynomial with known degree 
n, k +1 distinct roots, ro = 0,11,72,...,TR, with multiplicity mo,my1,me2,...Mx, 


respectively. Then the multiplicities are uniquely determined by co = 1,(,...,Ck- 


Furthermore, p may be determined by fewer than k proper coefficients when K is 


not algebraically closed. 


Theorem 5.3. Let p = Y.) a" € K [x] be a monic polynomial such that p(0) 4 0. 
Suppose p = [[j-, @" for q € K[az]. The multiplicity vector m = (m,...,m4)" is 
uniquely determined by the t proper leading coefficients if and only if V € K is 


non-singular where 


Vig= doom. 


rq; (r)=0 
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Remark 5.4. Observe that when gq; = x — 1; (i.e., p splits over K) Theorem 5.3 


provides the same conclusion as Theorem 5.2. 
We begin with a proof of Theorem 5.2. 


Proof of Theorem 5.2. Fix such a monic polynomial p with distinct roots rg = 0, 11, 


.., r~ With respective multiplicities mo, m1, ..., Mz. Ignoring ro for a moment, let 
r=(rj,...,7,) and m= (m,,...,m,)?. We denote the Vandermonde matrix 
1 il 1 
ry ue) Vk 
= 2 2 2 
4 ry "9 Vk 


and consider 


1 Lame Tk 
Pe eee r2 
1 19 k 
Vo = 
kok k 
Ty 9 ln 


Let p € K” where 
k 
Pis= Dorm, 
j=l 
then 


Vom = p. 


Notice that Vo = V diag(r) is non-singular as it is the product of two non-singular 


matrices. We have then 


m=V, 'p. (5.1) 


We present a formula for p which is a function of only the leading k + 1 coefficients 


of p. 
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Let A be the diagonal matrix where r; occurs m; times and note 
p(x) = det(al — A). 
By the Faddeev-LeVerrier algorithm (aka the Method of Faddeev [19]) we have for 
eel 
1 j ‘f Jj 
Cj = GF oy Chi tr(A’) = FG LB Cj-iPi- (5:2) 
i=1 i=1 


Leto = (1, C2, -..,¢)’, A = — diag(1, 1/2,...,1/k), and 


Co = 1 0 ... O 
Cy Co wink 

C= 
Ck—1 Chg .-. Cg 


By Equation 5.2, c = AC’p. Moreover, as A and C are invertible we have p = 


(AC)~‘c. It follows from Equation 5.1 that 


m = (ACV)~te. (5.3) 


Furthermore, mp = n — 1-m. 


We briefly remark about the proof of Theorem 5.2. Problem 5.1 has a flavor of 
polynomial interpolation: given k& points, how many (univariate) polynomials of de- 
gree n go through each of the k points? If n < k — 1 the polynomial is known to be 
unique and is relatively expensive to compute (as any standard text in numerical anal- 
ysis will attest). Our proof technique mimics this approach as the classical problem 


k-1 , .(k-1)-i 


of determining p = 37 cx , which resembles k distinct points {(2x;, y;)}*_,, 


can be solved by computing 


Vic=y 
where c = (Cg_1,Ck-2,---,C€0)', V is as previously defined given r; = x;, and y = 
(y1,--+;Yr)?. Suppose for a moment that each root is distinct so that 


k 


p(x) = [[(@ 1). 


i=1 


Al 


Then c,;, the codegree-j coefficient, is precisely the jth elementary symmetric polyno- 
mial in the variables r;,...,7rz. In the case of repeated roots we have that c; can be 
expressed using modified symmetric polynomials in the distinct roots where r; is re- 
placed with Gale . The expression for each coefficient via these modified symmetric 
polynomials is given by Equation 5.2. 

Note that if the roots of p are known, it is possible to determine p with fewer 
coefficients than the number of distinct roots (e.g., when p is a non-linear minimal 
polynomial). We modify Theorem 5.2 to include the case when some of the roots are 


known to occur with the same multiplicity. We now prove Theorem 5.3. 


Proof of Theorem 5.8. The proof follows similarly to that of Theorem 5.2. First 


suppose that V is non-singular. Let m = (m,,...,m,)7 and 
t 
pPi=>) De °) m; = (Vm); 
J=1 \r:q;(r)=0 


so that if V is non-singular, m = V~'p. Let A be defined as in Theorem 5.2: the 


diagonal matrix where the roots of q; occur m; times and note 
p(x) = det(aI — A). 

We have by the the Faddeev-LeVerrier algorithm, for 7 > 1 

i ee 

Cj =—-- os Cj-i tr(A’) =—-- S- Cj-iDi 

J j=1 J 421 

so that for c = (c1,...,c¢)7, A = —diag(1,1/2,...,1/#) we have 
m = (ACV) 'c. 


If instead m is uniquely determined by the first t proper coefficients then Vm = 


(AC)~'c has exactly one solution, hence V is non-singular. 


As a non-example, consider the minimal polynomial of a = V2, qa(x) = 27-2 € 


Z\z| and suppose 
p= qn = 24 as Or24-1 = Qdx24-? ea ey Z (2). 
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Observe that we cannot determine d given c; = 0; moreover, this conclusion is un- 


surprising, given the hypotheses from Theorem 5.3, since 
V = [v2 — v2] = [0] eC 


is singular. However, by inspection we could determine d given cz = —2d and in fact 
the matrix (/2" + (~V2)?] = [4] € C™? is non-singular. Indeed we could determine 


d with a simple change of variable: apply Theorem 5.3 to pp = (y—2)¢ where y = 2”. 


5.2. THE STABILITY OF COMPUTING MULTIPLICITIES 


We now consider the feasibility of computing m over C. In general, the matrix Vo 
in Theorem 5.2 may be poorly conditioned, so this calculation is often difficult to 
carry out even for modest values of k. The goal of this section is to show that if 
each root of a monic polynomial p(a) € Za] is approximated by a disk of radius at 
most €, a “reasonable” precision, then the interval approximating m,, resulting from 
a particular algorithm, contains exactly one integer. That is, we provide an algorithm 
for exactly computing m via SageMath [10] with substantially improved numerical 


stability over simply applying the Newton Identities. 


Theorem 5.5. Let p(x) = i_yca'' € Z{x] be a monic polynomial with distinct 


nonzero roots T1,...,T such that |ri| > |re] > --- > |rn| > 0. If each root is 


approximated by a disk of radius €¢ such that 


mr m n m n(1+o(1)) 
2 a 
© S pantIns (si) (Ga) 


where 
e M = max{max{|r; — 7;|}, 1} and m = min{min{|r; — r;| :7 A J}, 1} 
e R=max{|r,|,1} andr = min{|r,|, 1} 
e c= max{max{|c;|}",, 1}. 
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then the resulting disk approximating m; = (ACVo)~'c; € Z contains exactly one 


integer (i.e., the computation of m is exact). 


Notice that MRce > 1 and MRc = 2 for x” — 1 when n is even. Roots of 
unity occur frequently in the spectrum of hypergraphs; see Section 3. In particular, 
k-cylinders — essentially k-colorable k-graphs — have a spectrum which is invariant 
under multiplication by the kth roots of unity. Consider now p(x#) = 7” —1. We have 
m= V 2 — 2cos(2#) so that 
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(2 — 2.cos(2£)) "> 
D3n+7 5 


E< 


While ¢ may seem small, we are chiefly concerned with the number of bits of precision 
needed to approximate each root. Indeed for x” — 1 we need |Ine| = O(n Inn) bits 


of precision by the small-angle approximation. 


Remark 5.6. The bound on ¢ is “reasonable”, as the number of bits required to 
approximate each root is proportional to the number of distinct roots of p and the 
logarithms of the ratio of the smallest difference of the roots with the largest difference 


of roots, the largest root, and the largest coefficient. 


For the algorithm we describe below, the computational complexity of computing 
m is quasi-linear in nb where b = O(Ine) is the requisite number of bits of precision 
given. More precisely, the complexity is the maximum of O(nbIn nb) and O(n®) where 
6 is the best known exponent of matrix multiplication, currently [25] at 6 & 2.373 
but widely believed to be 2 + o(1). As we will see, m can be written as a product 
of n x n matrices and length n vectors populated by entries with O(nb) digits. The 
complexity follows as the number of digits needed to approximate each entry in the 
product of two such matrices/vectors is O(nb). 

In practice, the difficulty of computing m as described in Theorem 5.2 is in com- 


puting the inverse of the Vandermonde matrix, whose entries may vary widely in 
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magnitude and which may be very poorly conditioned. The task of inverting Van- 
dermonde matrices has been studied extensively. Eisenberg and Fedele [15] provide a 
brief history of the topic as well results concerning the accuracy and effectiveness of 
several known algorithms. However, these algorithms provide good approximations 
for the entries of V~', whereas we seek to express them exactly as elements of the 
field of algebraic complex numbers, since m is a vector of integers. Soto-Eguibar and 


Moya-Cessa [37] showed that V~' = AWL where A is the diagonal matrix 


n 1 tig ott 
Mates Peer aa ee 
0 OF jy 
W is the lower triangular matrix 


0 Eg 
Wig = 
Tkajtirei(Ti — Th) + otherwise, 
and L is the upper triangular matrix 
0 <4 
Lig =4 1 Spes) 


L; lj 1,-L; 1gjri-1 14 € [2,n],7 € [2,2 — 1]. 


Using this decomposition, it is possible to compute m exactly. To prove Theorem 5.5 
we first provide an upper bound for the diameter of the disk approximating an entry of 
A, W, and L, respectively; to do so, we extensively employ computations of Petkovi¢ 
and Petkovié [29] found in Chapter 1.3. We present the necessary background here. 
Let D(z,¢) be the open disk in the complex plane centered at z of radius ¢. For 


A = D(a,r,), B = D(b,r2) complex open disks, we have 
1At+tB=D(atb,r, +19) 
b r 
3. AB = D(ab, |alr2 + |b|r2 + rire) 
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In particular, for the special case of A” we have 
D(a,ri)" = D(a", (Ja] — ri)" — Jal"). 
Moreover, given 0 <r; <1 < |a| 
(Jal — ri)" — lal” < ri(2Ial)” 


since 


n n “ n n- n 
(Jal — ra)" — Jal” < > (i) tle * <ri(2lal)”. 
k=1 


Finally, let d(A) = 2r, denote the diameter of A and let 
Al = lal +r 
be the absolute value of A. Then for u € C we have 
1. d(A+ B) = d(A)+d(B) 
2. d(uA) = |uld(A) 


3. d(AB) < |Bld(A) + |Ald(B) 


(5.4) 


(5.5) 


For the remainder of this paper some numbers will be exact (e.g., rational num- 


bers) while others will be approximated by a disk. The non-exact entries of a matrix 


M € C”*?” will be referred to as disks; this will be clear from the problem formulation 


or derived from the computations. With a slight abuse of notation we use d(M;,;) 


and |M;,;| to denote the diameter and absolute value of the disk approximating the 


entry M;,. Moreover, we write 


d(M) = max{d(M,,;) : i,j € [n]} and |M| = max{|M,,| : 7,7 € [n]}. 


In the case when the entry is exact, the diameter is zero and the absolute value (of 


the disk) is simply the modulus. 
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Theorem 5.7. Assume the notations of Theorem 5.3, let V denote the Vandermonde 
matrix from the proof of Theorem 5.2, and let V-! = AWL by Soto-Eguibar and 
Moya-Cessa [37]. Then 


2n+4 n 
d(V~*) & 2 = (=) i 
m m 
and 
WV} <n (~—) | 
m 
Proof. Let 


dD; = Dri, é) 


denote the disk centered at r; with radius ¢. By Equation 5.5 we have for s 4 t 


aa) <4((5—) ) 


“ QE m se 
S? (= =: =) (= - =) 


gan+2 


ee 2 
— mnt2 53 


since e < m/4, 
d(W) <d((D,.—.D;)") <2" M"e, 
and 


d(L) < d(D") < (2R)"e. 


We first consider d(AW). Observe that AW is upper triangular and each nonzero 


entry of AW is a product of exactly one nonzero entry of A and W. In this way 


g2n+3 My” 
d((AW )ig) < |[Wigld(Ais) + [Avsld(Wig) < (—) E 


m2 \m 
and 
M n 
|AW| <2 (—) 
m 


We now determine d(AWL) by first computing 


gan+4 ey 
——— Ss 


m2 


d((AW)inLng) < |Legld(AWix) + |AWix|d(Lej) < 


m 


AT 


Hence 


- gantry RM te 
= 
d(V~1) = AWE) < max) d(AW) alas) $5" (>) § 
and 
ee < 2n (=*) 
m 


In our computations we are concerned with Vo = V - diag(r) where diag(r) = 
diag(ri1,...,7n) so that 
Vj. =diag(r)-'V—*: 


The following Corollary is immediate from the observation that 


Corollary 5.8. 


gQ2n+6 M n 
d(Vo}) <* (=) 
m2r m 
and 
2n (RM\” 
vc 28 (BMY 
r m 


We are now able to prove Theorem 5.5. 


Proof of Theorem 5.5. Recall m = V) 'C~!A~te as defined in the proof of Theorem 
5.2. Fortunately, the remainder of the computations are straightforward as C~!, A7?, 


and c have integer, and thus exact, entries. As 


1 = ee ree 
= ee Gg 2 


we have 


alV,: CG Sine” 


Pe aoe ai ee (MaeN". 


m2r m m2rc m 
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Further, since A~! = — diag(1,2,...,) we have 


g2n+6 4 MR n 
VC eave = ( “) é 


m2rc m 


and, finally, 
dV i Vela nemdVyanor ks) 


Oar oy ® (==) 1 
< E< 5" 


m?2r m 


Thus each interval will contain at most one integer as desired. 
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CHAPTER 6 


DETERMINING THE SET SPECTRUM OF A HYPERGRAPH 


In the previous chapter we showed that the characteristic polynomial of a hypergraph 
could be stably computed given its set spectrum and an appropriate number of leading 
coefficients. We now consider the problem of determining the set spectrum of H. In 
this chapter we summarize the Lu-Man Method which allows us to compute the 
totally nonzero eigenvalues (defined below) of a hypergraph. We then prove that 
the set spectrum of a hypergraph is contained in the union of the totally nonzero 
eigenvalues of all of its subgraphs. We further provide a partial characterization of 
precisely which of these eigenvalues extend to eigenvalues of the host hypergraph. 
This methodology is novel to spectral hypergraph theory as the same statement is 


not true for graphs (c.f. Cauchy Interlacing Theorem). 


6.1 THE LU-MAN METHOD 


The Lu-Man Method was introduced in [27] to study the spectral radius of a hy- 
pergraph and was further developed in [2], [40], etc. The method is based on the 
concept of an “a-consistent labeling” which is an assignment of complex values to 


every vertex-edge pair which satisfy certain conditions. 


Definition 6.1. ((40]) A generalized weighted incidence matrix B of a hypergraph H 
isa |V| x |E| matrix such that for any vertex v and any edge e, the entry B(v,e) = 0 


ifv ge and Biv,e) £0 ifu ce. 


Definition 6.2. ({40]) A hypergraph is called consistently generalized a-normal if 
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there exists a generalized weighted incidence matrix B satisfying 
1. Vevee B(v, e) = 1, for any v € V(A) 
2. Ince B(v, e) = @ for any e € E(H) 


3. for any cycle upeiv1e2... ve (ve = Vo), 


A totally nonzero eigenvalue of a hypergraph is a nonzero eigenvalue which cor- 
responds to an eigenvector with all nonzero entries. We denote the set of all totally 


nonzero eigenvalues of H as 0+ (H) C a(H). 


Lemma 6.3. (/40]) Let H be a connected k-uniform hypergraph and \ # 0. Then 
A is a totally nonzero eigenvalue of H if and only if H is consistently generalized 


a-normal with a = \~*. 


As the totally nonzero eigenvalues of a hypergraph are roots of its characteristic 
polynomial it is natural to define the a-polynomial of a hypergraph so that o+(H) = 
{A : ay(A) = 0}. When appropriate, we write a(H) for simplicity. We refer to the 
process of determining the a-polynomial of a hypergraph and all of its subgraphs via 
a-consistent labelings as the Lu-Man Method. 

Computing o(H) by way of ot(H) for H C H involves solving smaller multi- 
linear systems than the one involved in computing a(#). Generally speaking, |o(H)| 
is considerably smaller than the degree of ¢(H). In practice, this approach has 
yielded ¢(H) when other approaches of computing ¢(H) via the resultant have failed. 
In general, using the Lu-Man Method to determine the set of all totally nonzero 
eigenvalues of a hypergraph can be difficult as it involves determining a Grdbner 
basis a multi-linear system of (possibly non-homogeneous) equations with k|E| + 1 


unknowns; however, this method has been fruitful in practice. 
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We begin by providing a partial characterization of eigenvalues which are inherited 
from subhypergraphs. This characterization allows us to fully describe the spectrum 


of hypertrees. We conclude by providing a spectral characterization of power trees. 


6.1.1 INHERITING EIGENVALUES FROM SUBGRAPHS 


Recall that a vector is totally nonzero if each coordinate is nonzero and an eigenpair 
(A, x) is totally nonzero if X # 0 and x is a totally nonzero vector. Given a vector 
x € C” the support supp(x) is the set of all indices of nonzero coordinates of x. Let 
x° denote the totally nonzero projection (by restriction) of x onto C!**PP!, For ease 
of notation, we assume that the coordinate indices of vectors agree with the vertex 
labeling of the hypergraph under consideration. We denote the induced subgraph of 
H on U C V(H) by 


H|U] = (U, {v1... un € E(H) : v; € US) 


and write H C H to mean H = H[U] for some U C V(H). We begin by proving that 


an eigenpair of H can be projected onto its support. 


Lemma 6.4. Let (\,x) be a nonzero eigenpair of the normalized adjacency matrix of 


ak-uniform hypergraph H. Then (XA,x°) is a totally nonzero eigenpair of H|supp(x)]. 


Proof. As (\,x) is an eigenpair of H, 
n 
as HOed 
S- Qyiniz...ip VigXiz +++ Vi, = Ax; 


12.12,..,tRK=1 


for 7 € [n] by definition. Let m = |supp(x)| and suppose without loss of generality 
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that supp(x) = [m]. For 7 € [m] we have 


= S- Qjiziz...igVigliz .++ Viz 
ig.ig,...ip=1 

m 

S- Qjiziz...i,VigViz -- + Vi,- 
148s... 501 


Thus, (A, x°) is an eigenpair of H|m] by definition; moreover, (A, x°) is totally nonzero, 


as each coordinate of x° is nonzero by construction. 


It is not hard to see that if (A, x) is a nonzero eigenpair of H then \ is an eigenvalue 
of H|supp(x)| C H, the subgraph of H induced by the support of x [8]. This implies 
that the nonzero set spectrum of a hypergraph is contained in the union of the totally 
nonzero eigenvalues of all of its subgraphs (c.f. Cauchy’s Interlacing Theorem). That 


istor A 
aH) Go| se" (A). 
HCH 
We now consider the inverse problem: when can we extend a totally nonzero eigenpair 


of an induced subgraph of H to an eigenpair of H? 


Definition 6.5. An induced subgraph H|U] C H is isolated if E(H[U U {v}]) = 
E(H{U]) for all v € V(H). Moreover, a hypergraph is insular if every connected 


induced subgraph is isolated. 


Intuitively, an isolated induced subgraph has the property that including any one 


additional vertex to its vertex set does not induce additional edges. 


Lemma 6.6. Suppose H = H|U] C H is an isolated induced subgraph, then for 
(A,x) € ot (H) we have (A,x’) € o(H) where x!, = x, forv € V(A) and is zero 


otherwise. 
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Proof. Let H = H{[U]. Recall that (A, x) must satisfy 


k-1 _ 
i 


i= AL Uj Ujn °° Lj, 0 


{6,91,525-sJk—-1 tC H(A) 


fori € V(H). Let i € U and note that for {7, 71, Jo,..-,jx_1} € E(H) there is some 


t € [k — 1] such that 7, ¢ U which implies x;, = 0. Indeed 
eh Eg PON ap oe 0. 


Whence (A, x°) € ot (Hf), 


ae a 
pi = AtZ - > Dy gg Pe Dips 
{4,91,525-sJk-1}€E(H) 
ee al ee ree 
= Ar; — Dy Lj Lig Lip = 0. 


{4.91592 Jk-1}€E(A) 


Instead suppose i ¢ U. Note that any edge containing 7 is not an edge in H. Moreover, 
as H is isolated in H, every edge {i, j1,J2,.--,Jxr-1} ¢ E(H) has at least two vertices 


not in U. Indeed, there exists a t € [k — 1] such that x;, = 0. In particular, 


a a | 0 


so that p; = 0. 


Observe that Lemma 6.6 yields a sufficient condition for extending totally nonzero 
eigenvalues of subgraphs to eigenvalues of the ‘host’ hypergraph. We refer to totally 
nonzero eigenpairs which have this property as extendable in H. We have shown 
that all totally nonzero eigenvalues of an isolated subgraph are extendable; however, 
the converse of this statement is not always true. Consider the 3-uniform simplex, 
K. cm We claim that the simplex inherits eigenvalues from the single-edge despite the 
single-edge not being an isolated induced subgraph. In the case when a hypergraph 
is insular, its spectrum is exactly the union of the totally nonzero eigenvalues of its 
subgraphs. 

The previous argument can be applied to an arbitrary induced subgraph of H[U] 


if one can show that p; = 0 for i € U. This yields the following. 
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Corollary 6.7. A totally nonzero eigenpair of H{U] extends to an eigenpair of H if 


and only if p; =0 fori ¢ U as in the proof of Lemma 6.6. 


Remark 6.8. Recall that Theorem 3.27 expresses ¢a(H) as a weighted sum of graph 
statistics. It is natural to question if our approach to spectral hypergraph theory, via 
hypermatrices and polynomial maps as discussed in the Introduction, is equivalent 
to that of a linear spectral theory in some way. Corollary 6.7 suggests that these 


two approaches are inherently different as demonstrated by this spectral inheritance 


property. 


In the following section we characterize the totally nonzero eigenvalues of a hy- 


pertree. 


6.2. SPECTRA OF HYPERTREES 


The following beautiful result was shown in [40]: the set of roots of a certain matching 
polynomial of a k-uniform hypertree (an acyclic k-uniform hypergraph) is a subset of 


its homogeneous adjacency spectrum. 


Theorem 6.9. (//0/) Let H be a k-tree. Then 


o(H) = S(-1)'|Miln 


i=0 


where M, is the collection of all i-matchings of H. 


We present the a-polynomials of 3-trees with six or fewer edges in Table 6.1 and 
the number of 3-trees is given by A003081 [36]. We now show how to obtain all of the 
eigenvalues of a hypertree, and use this description to give a spectral characterization 
of “power” hypertrees (defined below). We extend Theorem 6.9, as follows, to describe 


the spectrum of a hypertree, answering the main open question in [40]. 
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Lemma 6.10. A 3-graph H is insular if and only if it is a hypertree (i.e., H is 
connected and every nontrivial walk vie, ...en—1v1 must have a repeated edge). In 


particular, 


o(H) = {0}U LU o* (A) 
HCH 


where the union is taken over all connected induced subgraphs, if and only if H is a 


hypertree. 


Proof. Let H be a hypertree and let U C V(H) such that H = HU] is a subtree of 
H. Suppose e = {u,v,w} € E(H) for u,v € U and w € U’. Note that a hypertree is 
linear (i.e., two edges intersect in at most one vertex). It follows that e is the only 
edge of H which contains both u and v. Since u,v € U and H is connected, there 
exists a nontrivial walk from u, €1, v2,...,€:-1, v in H where all the edges are distinct. 
Note that wu, e1,V2,...,e€¢-1,U,e€,u is a cycle in H, a contradiction. Indeed no such 
e € E(H) exists so that H is isolated in H. 

Let H be an insular 3-graph. Note that H is linear, for if this were not the case the 
hypergraph would contain an edge which was not isolated. Moreover, H is acyclic as 


a linear 3-uniform hypercycle (a.k.a. a loose 3-cycle) is not insular as the path formed 


by removing one vertex is not isolated. 


6.2.1 SPECTRA OF POWER TREES 
The following generalizes the definition of powers of a hypergraph from [24]. 


Definition 6.11. Let H be an r-graph for r > 2. For any k > r, the k'" power of 


G, denoted H*, is a k-uniform hypergraph with edge set 


E(H*) ={eu {sixes Ve Kory SOE (GI, 
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Table 6.1. The a-polynomial of all 3-trees with six or fewer edges. 


E(T) a(T) 
[1,233] ge —1 
[E9- 3), 104.5) oD 
[152 3); 4,5) 6F 7] qo = 3 
[i2: 31 A, an [2, 6, 7] = Be 1 
[1, 2, 3], (1, 4, 5], (1, 6, 7], [1, 8, 9] o—A4 
erie San 8,9], 2, 6,7] x® — 4r3 + 2 
[1, 2, 3], [1, 4, 5], [2, 6, 7], [3, 8, 9] x? — 4x® + 3x° — 1 
[1, 2, 3], [1, 4, 5], [2, 6, 7], [4, 8, 9] gos Ag is3 
1, 2,3], (1,4, 5], (1, 6, 7], [1, 8, 9], [1, 10, 11 p35 
POs ae: 1.8.9], 1 10511) 1256.77 gp — 52°43 
1, 2,3], (1,4, 5], (1, 8, 9], [2, 6, 7], [2, 10, 11 p= be tA 
1, 2, 3], [1,4, 5], [1, 8, 9], [2, 6, 7], [3, 10, 11 go = 5g 5g aD 
1, 2, 3], (1,4, 5], [1, 8, 9], [2, 6, 7], [4, 10, 11 eo ob Loe 1 
1, 2, 3], (1, 4, 5], [2, 6, 7], (2, 10, 11], [4, 8, 9 x° — 5a +5 
1, 2, 3], (1,4, 5], [2, 6, 7], [3, 10, 11], [4,8,9 ge bg? 659 
1, 2, 3], [1, 4, 5], [2, 6, 7], [4, 8, 9], [6, 10, 11 Pen ae a oe ee 
19 3c, Sl y(1, 6; 7s (1: 8, 0); [et], hs 1913 x —6 
1, 2, 3], [1, 4, 5], [1, 8, 9], [1, 10, 11], [1, 12, i [2, 6,7 x? — 62° +4 
1, 2, 3], (1, 4, 5], [1, 8, 9], [1, 10, 11], [2, 6, 7], [2, 12, 13 x° — 62? +6 
1, 2, 3], [1, 4, 5], [1, 8, 9], [1, 10, 11], [2, 6, 7], [4, 12, 13 BP She eit =! 
19:3) (1, As5)s [08,9], 11, 12,1312; ao x? — 62° + 7x3 — 3 
1, 2,3], {1, 4, 5], [1, 8, 9], [2, 6, 7], ia 3,10, 11 x? — 6x® + 823 — 4 
13 2,3) 11; 4:51, [1s 8201; 122627). 12.1% 13), 10, 11. x? — 62° + 8x3 — 2 
1, 2, 3], [1, 4, 5], [1, 8, 9], [2, 6, 7], [3, 12, 13], (4,10, 11] | 21? — 6x9 + 97° — 52741 
1, 2, 3], [1, 4, 5], [1, 8, 9], [2, 6, 7], [4, 10, 11], (6, 12, 13 go be? 07" 38 
1, 2,3), {1, 4, 5], [1, 8, 9], [2, 6, 7], [4, 10, 11], (8, 12, 13 x? — 62° + 923 — 4 
1, 2,3], {1, 4, 5], [1, 8, 9], [1, 12, 13], [2, 6, 7], (6, 10, 11 o — br? +7 
1, 2,3], {1, 4, 5], [1, 8, 9], [2, 6, 7], (6, 10, 11], (6, 12, 13 x® — 623 + 8 
1,233), [045 5), (2,6, 712.12; 13), (3,10; 11), (4, 829 x? — 62° + 923 — 4 
1,953), 4,5), (2, 6.7 218,10 11) (428-0) 8 x? — 62° + 9x3 — 3 
L238) {EAs 5ly 256.719 13.10 11 (48 Ol Se 12 8h | et = Cae? 10a be? 1 
1,93), 114s Sls (D677 o (3:10. 114, 8-0 16,19: 13 x? — 6x® + 102? —5 
1, 2,3), {1, 4, 5], [2, 6, 7], [3, 10, 11], [4, 8, 9], (8, 12,13] | «!? — 6a® + 102° — 523 +1 
1, 2, 3], [1, 4, 5], [2, 6, 7], [4, 8, 9], [4, 12, 13], [6, 10, 11 x? — 62° + 925 — 2 
1, 2, 3], [1, 4, 5], [2, 6, 7], [4, 8, 9], [6, 10, 11], [8, 12, 13 x? — 6x® + 102? — 4 


57 


and vertex set 
V(H*) = V(G) = V(G)U {i.e € E(G), 9 € [k—-r]}. 


In other words, one adds exactly enough new vertices (each of degree 1) to each 
edge of H so that H* is k-uniform. Note that, if k = r, then H* = H. Adhering 
to this nomenclature we refer to a power of a 2-tree simply as a power tree. In this 


section we prove the following characterization of power trees. 


Theorem 6.12. Let H be a k-tree and let o(H) denote its (multiset) spectrum. Then 


o(H) C R[G] if and only if H is a power tree, where &, is a principal k™ root of unity. 
We recall the following Theorem from Cooper and Dutle. 


Theorem 6.13. /12/ The (multiset) spectrum of a k-cylinder is invariant under mul- 


tiplication by any k" root of unity. 


One can show by straightforward induction that a k-tree is a k-cylinder, so its 
spectrum is symmetric in the above sense. The following result, from [41], shows that 


power trees have spectra which satisfy a much more stringent condition: they are 


cyclotomic, in the sense that they belong to R|¢]. 


Theorem 6.14. /41/ If \ £0 is an eigenvalue of any subgraph of G, then \?/* is an 


eigenvalue of G® for k > 4. 


We restate Theorem 6.14 with the additional assumption that the underlying 
graph is a tree; the proof is easily obtained by applying Lemma 6.10 to the proof of 


Theorem 6.14 appearing in [41]. 


Corollary 6.15. If \ 40 is an eigenvalue of any subgraph of a tree T, then d?/* is 


an eigenvalue of T*® for k > 3. 


Note that Theorem 6.12 provides a converse to Corollary 6.15 in the case of power 


trees. In particular, appealing to Lemma 6.10, Corollary 6.15, and the fact that the 
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spectrum of a graph is real-valued, we have that the spectrum of a power tree is a 


subset of R[¢;,]. All that remains to be shown is that if a k-tree is not a power tree 


then it has a root in C \ R|¢;|. To that end, we introduce the k-comb. 


Let Comb, be the k-graph where 
Comb, = ((k*], {[k] U{{it+ th: 0<t<k—1}:é€e [k]}}}. 


We refer to Comb, as the k-comb. By the definition of power tree, a non-power 
tree H must contain an edge e incident to a family F consisting of at least three other 
edges which are mutually disjoint. This edge e, together with F, form a connected 
induced subgraph H' of H which is the k power of a t-comb for t = |F| > 3. It is 


straightforward to see that 
an (©) = Agomps() = comb, (x**), 


since matchings in H’ are simply k* powers of matchings in Comb;,; therefore, roots 
of a(H’) are k™ roots of reals if and only if the roots of a(Comb;) are ¢‘" roots of 


reals. We presently show that the spectrum of the k-comb is not contained within 


the k'* cyclotomic extension of R, completing the proof of Theorem 6.12. 


Lemma 6.16. There exists a root of a(Comb,) for k > 3 such that A € C \ RIG]. 


Proof. Let H be a k-comb where k > 3. By a simple counting argument, 
i=0 
where a = x*. Appealing to the binomial theorem we have 
a(H) = (1—a)* —a*?, 
Let 8 =a. Setting a(H) = 0 yields 


(6 -1)* = 8. (6.1) 


It is easy to see that (6.1) has precisely one solution when k > 3 is odd and precisely 


two solutions when it is even. In either case, as the number of solutions is strictly less 


than k it follows that there must be a non-real solution and the claim follows. 
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CHAPTER 7 


COMPUTATIONAL NOTES AND EXAMPLES 


We consider the problem of computing the characteristic polynomial of a uniform 
hypergraph. For a k-graph on n vertices the characteristic polynomial is the resultant 
of n homogeneous multi-linear polynomials of degree k —1. While there are numerous 
methods for solving such a system, computers must be employed even for the simplest 
of cases. In practice, built-in methods for solving such systems, via Sage [10] and 
Macaulay2 [21] for example, have proven to be insufficient. This motivates the need 


for a new approach which we now present by invoking our aforementioned results. 


1. Determine S, the set of induced subgraphs of H up to isomorphism. 


2. We say that H € S is a spectral candidate in H if there exists an H'’ C H 
such that H’ = H and H’ satisfies Corollary 6.7 (this includes the case when 
it cannot be determined if H’ satisfies the Corollary). Let S’ C S' be the set of 


spectral candidates of H. 


3. Compute a(H) for H € S" by way of the Lu-Man Method. The authors use Sage 
[10] to create the system of equations and solve the system using Macaulay2 
[21] by computing the generators of the Grébner basis of the ideal generated by 


the aforementioned equations. 


4. Let {p;})_, be the collection of minimal polynomials of a(H) for H € S’. Note 
that 2-"°¢(H) = [Ii_, p”. Furthermore, it may be the case that m; = 0 if p; 


is a factor of a(H) for H € S’ where Corollary 6.7 could not be verified. 
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5. In the case when o(H) is symmetric by some principle root of unity make an 
appropriate change of variable to simplify further computations. This symmetry 


will be clear from the minimal polynomials. 


6. Choose to apply Theorem 5.2 or Theorem 5.3. Generally speaking, Theorem 5.3 
is more efficient in the case when there are few minimal polynomials of larger 


degree. 


7. Determine an appropriate number of leading coefficients of ¢(H) (with the 
change of variable, if applicable) using Theorem 3.27. Note that the coefficients 
can be computed more efficiently by taking the structure of the hypergraph into 


account in order to reduce the search space of the Euler rootings. 


8. Apply Theorem 5.2 (or 5.3) to compute the multiplicities of each root (or min- 


imal polynomial, respectively). 


We now demonstrate this algorithm by computing the characteristic polynomial 


of the Hummingbird hypergraph and the Rowling hypergraph. 


7.1 THE HUMMINGBIRD HYPERGRAPH 
Consider the hummingbird hypergraph B = ((13], £) where 
Beat F253 }34 4, oh (ly 6s Ch 425 8, 9} 3) 10, dS To Tah 


We present a drawing of 6 in Figure 7.1 where the edges are drawn as shaded in 


triangles. Note that 
deg(#(B)) = n(k — 1)""' = 13 - 24 = 53248 
and, since B is a hypertree (and thus a 3-cylinder), its spectrum is invariant under 


multiplication by any third root of unity [12]. Appealing to Lemma 6.10 we have that 
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Table 7.1 Subgraphs of 6 and their a-polynomial. 


Tep a(T) 
P= 8; ge 41 
Py = Sp 79 

53 eas 


(10.8) Abe: BV ACT 67412. 8.9 4) ro — 473 +2 


{{1, 2,3}, {1, 4,5}, {2, 8, 9}, {3, 10,11} | 2? —42°+ 32° -1 
B— {2,8, 9} g® — bar? +4 

B — {3, 12, 13} php £5 = 9 

B x? — 62° + 87? — 4 


S’ is the set of all subtrees of 6, up to isomorphism. We compute the a-polynomials 


of subtrees of B via the Lu-Man Method. 


6(B) =a" (x° — 62° + 82° — 4)™ (2° — 52° + 5a? — 2)™ 
(a? = 1)" (2° — 42? + 2)™ (2? =—4z® + 32? = 1) 
- (2° — 3a? + 1) (2? — 3)" (x? — 2). 


With the intent of applying Theorem 5.2 to ¢(B) we consider the change of variable 


y = x? and observe that we need to determine c3,¢¢,...,C4g as there are sixteen 


63 


distinct nonzero cube roots. Appealing to Theorem 3.27 we compute 


cz = —18432 

ce = 169843968 

cg = — 1043209971456 

C12 = 4804960103034624 

C15 = —17702435302276375440 

Cig = 54341319772238850901668 

C21 = —142960393819753656566577552 

Coq = 329036832924106136747171871042 

C97 = —673063350744784559041302787109576 

C39 = 1238925078774563882036470496247467682 

C33 = —207289173587094969593028654258099 1559916 

C36 = 3178738418917825954994865036362341584776658 

C39 = —4498896549573513724009044022281523093964642496 
C42 = 5911636016042739623328802656744094043553245557890 
C45 = —7249053168113446546908444934275696322928768819713512 


Cag = 83322309372 1367842625849 1158832963153453272812465162851 


Using Theorem 5.5 we have M < 3,m > .14,R < 4.39,r > .38, and c = cag so that 
each root of ¢(B) needs to be approximated to at most 3091 bits of precision. Using 


SageMath [10], we obtain 


gp = 2783 (79 — 6a + 8x3 — 4)9(49 — 52% + 5a — 2)97 
- (x8 — 1)1782(8 — 4g + 2)488 (29 — 46 4 303 — 1)8%4 


; (x° = 3,3 Bs Lye? = Neal eos = Qe 
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In Figure 7.1 we provide a plot of o(B) drawn in the complex plane where a disk is 
centered at each root and each disk’s area is proportional to the algebraic multiplicity 


of the underlying root in (B). 


Figure 7.1 The hummingbird hypergraph and its spectrum. 


7.2. THE FANO PLANE AND ITS SUBGRAPHS 


We compute the characteristic polynomial of the Rowling Hypergraph, that is, the 
Fano Plane with two edges removed. The characteristic polynomial of the Fano 
Plane and the Fano Plane with one edge removed remain to be determined. In this 
chapter we present our progress towards computing the characteristic polynomial of 
the Fano Plane with one edge removed and then the Fano Plane with the intent of 
demonstrating the current computational limits of this method. We encourage the 
eager reader to improve upon them. Before we proceed, we present the first few 


leading coefficients of the aforementioned hypergraphs in Table 7.2. 
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Table 7.2 Leading coefficients of the Fano plane and its subgraphs. 


FP-2 FP-1 FP 
Co 1 1 1 
Cr 0 0 0 
C2 0 0 0 
C3 -240 -288 -336 
C4 0 0 0 
C5 0 0 0 
C6 28320 40788 55524 
C7 0 0 -696 
Cg 0 0 0 
C9 -2190860 -3788016 -6017746 
C10 0 0 220038 
C11 0 0 0 
C12 125012034 259553826 481293561 
C13 0 0 -34237560 
C14 0 0 -122004 
Ci5 | -0612445168 | -13997317932 | -30303162330 
C16 0 0 ? 
Ci7 0 0 ? 
Cig | 206518037092 | 618904026740 i 


7.2.1 


Consider the Rowling hypergraph,! 


R= UA 3 14 op 6 Th 2oh Gh 3s5, CP: 


Note that R is the Fano Plane with two edges removed. A drawing of R is given 
in Figure 7.2 where the edges are drawn as arcs and its multi-spectrum is drawn 


such that each disk is centered on an eigenvalue of R and the area of each disk is 


THE FANO PLANE WITH Two EDGES REMOVED 


proportional to the multiplicity of the approximated root. We have 


deg(¢(R)) = n(k —1)""' =7- 2° = 448. 


'The name was chosen for its resemblance to an important narrative device[32]. 
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There are four induced subgraphs of R, up to isomorphism, and we pick a represen- 


tative of each class so that S = {H;}4,: 
A, = RI [3] = P,, Ao = R5] = P», H3 = R{|7| \ 4], Ha = R. 


Again we have that H, and Hy, are both isolated subgraphs of R so that by Lemma 
6.6 their totally nonzero eigenvalues extend to 7. 

We claim that the totally nonzero eigenvalues of Hj and H3 are not extendable 
in R. Adhering to the language of Lemma 6.6 we set U = [5]. Consider R[U U {6}| 
and observe 

Pe = ANG — oes = —Loits 
since 7g = 0 as 6 ¢ U. However, because 2,5 € U we have that pg 4 0. By Corollary 
6.7 the eigenvalues of H» are not extendable to R. This argument can be generalized 


to prove the following. 


Lemma 7.1. [fH is ak-graph and H|U] has the property that there exists v ¢ U such 
that |E(H|U U v})| — |E(H[U])| = 1 then the eigenvalues of H|U] are not extendable 
in H. 

Lemma 7.1 implies that the totally nonzero eigenvalues of H3 are not extendable 
in R with the choice of v = 4. 

Indeed, S’ = {H,, H,} so that o(R) = ot(P,) Uot(R). We have previously 
shown that a(P;) = x? — 1 and we have by the Lu-Man Method that 


a(R) = (a — 132 + 650° — 147° + 1572? — 64) (2° — 2 + 2) (2° — 172° + 64). 


Notice that R is not a 3-cylinder; however, its spectrum, like that of 3-cylinders [12], 
is invariant under multiplication by any third root of unity (see Lemma 3.11 of Fan, 


et al. [17]). We have then that 
WR) =e" (2? = 1)™ (2 — 132" + 652" — 1472" + 1572° — 64)" 


(gi =a? 42) (2° = 172? +64)", 
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With the intent of applying Theorem 5.3 to the minimal polynomials of 
x MR ba (all*), 


we verify that 


1 39 -3 161 
1 103 -—5 1649 
1 ‘8 1 T7229 
is non-singular. Indeed, we need to determine the first four proper leading coefficients 


of 2~°/3dp (2/3), or equivalently c3,c¢, C9, C12 of dx. We have 


C3 = —240 
Ce = 28320 
Cg = —2190860 


By Theorem 5.5 we have M < 4.5,m > .69, R < 2.25, and r = 1 so that at most 252 


digits of precision are required to approximate each root. We compute 


bp =x (a? = 1)" (e? = 182"? + 652° — 14728 £ 15727 — 64)" 


(a= 2? +2)%(2° = 172? + 64)? 


7.2.2. THE FANO PLANE WITH ONE EDGE REMOVED 


Let 
FS (754 $152,34541,6, 6154258 Ppt 2a 0 a has oy Ch: 


be the Fano Plane with one edge removed. There are four spectral candidates we 


need to consider: 
Ay =F ||3\| = Py fo =F fl, 25955, }| = Po sg [C4 a ae 
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Figure 7.2. The Rowling hypergraph and its spectrum. 


Again we have that H, and Hy, are isolated subgraphs so that their totally nonzero 
eigenvalues extend to F. Observe that Hz and H3 are not isolated subgraphs and 
further the assumptions of Lemma 7.1 are not satisfied. Thus, we cannot conclude 
whether or not the totally nonzero eigenvalues of Hj or H3 extend to F’. In this case, 
we assume Hy and H3 are spectral candidates of F' so that S’ = S. By the Lu-Man 
Method we have 


a(H,) =z —1 
a(H) =z —2 
a(H3) =(a? — 8)(x® — 52° + 8) 
o(H4) =(x? — 2)(x° — 18)(2° + 2? +1) 
(x® — 172° + 80)(x" — 62° + 192° — 192? + 8). 
Observe that a(H2) and a(H4) have a common root. We have then that 
br(a) =a? (2® = 1)™ (2? = 2)" (2? = 8)" (2® — 52? 48)™ 
~ (g? = 18) (2° +a? + 1) (2° = 172° +80)" 
-(¢? — 62° + 197° — 1927 + 8). 


In order to determine mo,...,™mg we would apply Theorem 5.3 which requires 
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knowing the first eight proper leading coefficients of 2~"°/3¢p(x!/3). As shown above, 
we know only the first six proper leading coefficients. Unfortunately, our current 


routines are not able to compute the remaining two coefficients due to run time. 


7.2.3. THE FANO PLANE 


Consider the finite projective plane of order 2, or more affectionately, the Fano Plane 
F = (ttl 2.3414. ob 41.6, 42/4 42.5, 6h 3, 406) 43.5, Th): 
Once again, there are four spectral candidates of F, 
y= F|3|| Pits =F |[5|| ea Fl Va ae 


In this case, we have been unable to compute a(F). We have attempted to use several 
different built-in methods in Macaulay2, all of which failed due to run time. Once 
one determines a(F) the task of determining the multiplicities remains. It is unclear 
if we will be able to apply Theorem 5.3 or if we must use Theorem 5.2. If we apply 
Theorem 5.2, the first fifteen coefficients are required before considering a(F)! This 
means we will need to compute a number of coefficients equal to the degree of a(F). 
It appears that considerable work must be done in this direction. Not just in terms 
of being able to compute a sufficient number of coefficients, but also in computing 


the a-polynomial of an arbitrary hypergraph. 
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CHAPTER 8 


OPEN PROBLEMS 


We conclude by presenting open problems related to our work. 


8.1 COSPECTRAL HYPERGRAPHS 


Two graphs are cospectral, or isospectral, if they have the same multi-set spectrum 
(i.e., they have the same characteristic polynomial). In 1973, Schwenk showed that 
almost all trees have a cospectral mate [34]. Thirty years later, van Dam and Haemers 
published a survey paper on graphs which are uniquely determined by their spectrum, 
abbreviated DS [13]. In it, they suggest that it is conceivable that almost all graphs 
are DS. In [4], Bu, Zhou, and Wei show that complete k-uniform hypergraphs are DS 
(so the complete graph is DS for any uniformity). Determining which hypergraphs 
are DS appears to be an interesting question and it is unclear a priori if this question 
has a different fate than the graph case. 

To demonstrate this, recall that the smallest pair of cospectral simple graphs is the 
star on five vertices and the 4-cycle with an isolated vertex. We show that the power 
graphs of these two graphs are no longer cospectral. Let S and C be the 3-graphs 
formed by adding a unique vertex to each edge in Ss and Cy U {v}, respectively. We 
will show that S and C are not cospectral despite their base graphs being so. From 


the Lu-Man Method we have 


HS) =a (a = 1)™ (a — 2) (a" = 3) (a —4)™ 
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where m; > 0 and further 
H(C) = "(29 — 1)*(0° — 2)"%(a9 — 4 


where n; > 0. Observe that a(S) and o(C) are not equal as sets. Indeed S and C are 


not cospectral mates. This begs the question: 
Question 8.1. What are the smallest cospectral k-graphs for k > 3? 


In order to answer this question, one would first determine two hypergraphs which 
have the same set spectrum and then one would compute the characteristic polyno- 
mial of each hypergraph, which is intractable in general. This motivates the following 
variants of cospectral hypergraphs. Two hypergraphs are weakly cospectral if their 
spectra are equal as sets. Moreover, if H, and H2 are hypergraphs we say that Hj is 
subspectral to Hy» if ¢(H,) | o(H2). Under these definitions we can restate a conjecture 


of [8] as follows. 


Conjecture 8.2. A hypertree is subspectral to any hypertree which contains it as a 


subtree. 


We conclude this section by presenting an example supporting this conjecture. 
Consider the 3-uniform hypergraphs 
Fy = (191; 4 1,.2:3)1L4, Th 2,-558} 413,69} }) = Combs 
H2 = ({11], gee 2, 3}, {1, 4, ide {3, 6, OF; (i; 10, 11}}) 
Hs = ((11], EH1) U E(H2)). 


Figure 8.1 gives a drawing of H3 (the striped subgraph is H;) and a plot of the roots 


of ¢(H3), with a circle centered at each root in the complex plane whose area is 


(2 


proportional to the multipli 
o(H1) = De Cod = 4x® 7 
(Hz) = Dad Cd a 4x? 7 


(Hs) = 23729 — 56 


: (x° -_ 3,3 ae 1 ee? = Se = ae oo = 


city of the root. We have computed 


L 3,3 _ 1a? = 373 at Lye? = Qe? _ Mee 


L 2) (2° a 37° As Lyte? = 3)?" (x3 = Dd eo = De 
a 53 = 2S (x? 2! Ar® a 373 fe! 1) Ge = Ar? a gy 


iD ae 


Note that H, and Hy» are subspectral to H3. 


Figure 8.1 Hz3 and its spectrum. 


8.2 MULTIPLICITY OF THE ZERO EIGENVALUE 


A characterization of the multiplicity of the zero eigenvalue for the adjacency char- 
acteristic polynomial of a graph remains open. For convenience let mo denote the 
multiplicity of the zero eigenvalue for a given polynomial. Notice that one can provide 
an upper bound on mo by showing that a particular coefficient of ¢(H) is non-zero. 
Combining this idea with the Harary-Sachs Theorem gives the only known result in 
this direction for the adjace 


uniform) tree then mp is the size of the largest matching of T. Note that the same 
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ney characteristic polynomial of a graph: if T is a (2- 


statement for hypertrees cannot be true (consider Lemma 6.10); however, we know 
that the spectrum of a hypertree can be written in terms of the matching polyno- 
mials of all of its subtrees. Indeed, if one could characterize the multiplicities of the 
matching polynomials of a given subtree of a hypertree one could answer the following 


question. 


Question 8.3. Express the multiplicity of the zero eigenvalue of a hypertree in terms 


of the sizes of its matchings over all of its subtrees. 


In our proof of Lemma 3.29 we showed that, for 2-graphs, the summands in ¢4(G) 
arising from Veblen graphs which are not elementary graphs necessarily summed to 
zero. We define the coefficient threshold of ¢(H) as the least codegree at which the 


coefficients of ¢(H) cancel thusly. 


Definition 8.4. For an integer v > 0, the coefficient v-threshold of a k-graph H, 


denoted Th,(H), is the least integer such that for d > Th,(H) 


S>  (-((k-1)"))Cu(#H CH) =0. 


Hevi(H) 
Note that the contribution of H to the codegree-d coefficient of (G) for H C G is 
zero if d > Th,(H) where v = |V(G)|. As an example, we show that the v-threshold 


of the 3-uniform edge is 9 - 2°~°. 
Lemma 8.5. Let e be the 3-uniform edge. Then Th,(e) =9-2°~° for v > 3. 


Proof. For n > 0 define f,(t) : Z+ + Q by f,(0) = 1 and 


frlt)= YO (-(2"))Cu(#H C e),t > 0. 


HEV3,(e) 
Observe that f3(t) = @3:(e) by Theorem 3.27. We conclude by showing f,(t) = 


(-1)'(?7""). Considering the characteristic polynomial of a single 3-uniform edge, 
f3(0) =1, fa(1) = —3, fs(2) =3, fa(3) =—1, fa(4) =0 
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and f3(t) = 0 for t > 4. Indeed f3(t) = (-1)'(3) for all t. Suppose that for all n, up 
to some fixed n, we have f,(t) = (-1)'(?7"") for all t. Consider f,41(t). We claim 


t 


Fnoi(t) = 2s Pit =F): 


Recall that the associated coefficient is multiplicative over components. It follows 


that 
t 
2 In ID Fn( Ce j) — SS yo (—(2")) Cp, ' (—(2"))) Cx, 
j=0 HEV3(e) H=H,UH2 
Hi NH2=0 
7 S- oe 2) ye NDOy = Frsi(t). 
HEV 3 (e 
We have 


Ant EOE MET) 


where the first equality follows from the inductive hypothesis and the second equality 
is given by the Chu-Vandermonde identity [6]. As Th,(e) = f.(t) we have that 


‘Thyte) = 0-00" asf (3-29) = tort = 3-2" and jf. = 0 for f= 372". 


Conjecture 8.6. Jf H C G are k-graphs, with k > 2 where |\V(G)| = n then 
Lh) <7 Pin Gi 


This conjecture implies 
mo < deg(¢(G)) — Th, (H). 


Note that the restriction of k > 2 is necessary as the conjecture is not true for graphs. 
For example, Cg C Ke¢ and Th(K6¢) = 2 < Th(Cs) = 6. More generally, it would 
help our understanding of the multiplicity of 0 to have a better understanding of 


Th,(H), and so we ask: 


Question 8.7. Show how to compute or estimate Th,(H) for various hypergraphs 
Fh 
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8.3 COMPUTING THE Q-POLYNOMIAL OF A HYPERGRAPH 


The Lu-Man Method is quintessential to our algorithm for computing the character- 
istic polynomial of a hypergraph. For this reason it is of great interest to the authors 
to be able to compute the a-polynomial of a given hypergraph. At the present, the 
authors have had mixed success in this direction. Recall that computing a(H) in- 
volves solving a system of equations with k|E(H)| + 1 unknowns with a number of 
equations equal to the sum of the sizes of the vertex set, edge set, and cycle basis. 
In general, we have found that the a-polynomial of a hypergraph with a vertex (or 
vertices) of degree 1 can be computed efficiently. For example, the a-polynomial of 
a tree is its matching polynomial, but the a-polynomial of the Fano Plane remains 


unknown. 
Question 8.8. Show how to compute the a-polynomial of a hypergraph efficiently. 


Another approach would be to characterize the a-polynomial of a certain family of 
hypergraphs. A natural family to consider would be that of loose-cycles. We present 
the a-polynomials of loose 3-cycles with 12 or fewer edges in Table 8.1. Notice that, 


for the values given, a(C) | a(c?) given i | j. 


Table 8.1 The a-polynomial of loose 3-cycles with i edges where y = 2°. 


i aC!) 

3 (y—1)(y—4) 

4 (y —2)(y — 4) 

5 (y — 4)(y* — 3y + 1) 

6 (y —1)(y— 3)(y—4) 

7 (y — 4)(y* — 5y? + 6y — 1) 

8 (y — 2)(y — 4)(y* — 4y + 2) 

9 (y — 1)(y — 4)(y? — 6y* + 9y — 1) 

10 (y — 4)(y* — 5y + 5)(y* — 38y + 1) 

11 | (y — 4)(y> — 9y? + 28y° — 35y? + 15y — 1) 
12 | (y—Iy—-2)(y—-3)y¥—-4(y* — 4y +1) 
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Question 8.9. Characterize the a-polynomial of loose 3-cycles. 


Aside from computing the a-polynomial of a hypergraph we would like a better 


test for determining if H|U] C H is a spectral candidate of H. In particular, in the 


case when |E(H)| — |E(H[U U {v}])| > 1. Such a test would permit the study of hy- 
pergraphs with few eigenvalues. Famously, strongly regular graphs are characterized 
by having few distinct eigenvalues. A natural question to ask is whether the power of 
a strongly regular graph retains this property and an answer to this question would 


serve as a basis for answering the following. 


Question 8.10. Characterize k-graphs which have few distinct eigenvalues. 


8.4 ‘THE SPECTRA OF A RANDOM HYPERGRAPH 


A central topic in spectral graph theory is the spectra of random graphs. This 
area addresses questions concerning quasirandomness, graph expansion, and mixing 
time of Markov chains. The spectrum of an (Erdés-Rényi) random graph is well 
understood; however, understanding the spectrum of a random hypergraph is a more 
subtle matter. We present two approaches to understanding the spectrum of a random 
hypergraph. We first address the problem given an analogue of the Weyl inequality, 
then we address the problem from the perspective of quasirandom hypergraphs. 
Consider for a moment the simpler boundary case of this question where each 
edge is chosen with probability 1. Here, the resulting hypergraph is the complete 
k-graph on n vertices, denoted kK). We could approximate the spectrum of K‘) by 
comparing it to a suitable hypermatrix if we had an analogue of the Weyl inequality 
(as stated below) for symmetric hypermatrices. Consider the normalized all-ones 
order n dimension k hypermatrix, (k — 1)!~!J* and observe that this hypermatrix 
agrees with A(K“)) except for tuples of indices whose coordinates are not all distinct. 


In [11], Cooper showed that there is an approximate bijection between the set of 


ws 


eigenvalues of K‘”) and the set of eigenvalues of (k — 1)!-'J* for the case of k = 2,3. 
More precisely, the set of eigenvalues L of K() and the set of eigenvalues M of 
(k —1)!"'J* satisfy 6(L,M) = o(n*~') where 6 is the Minkowski distance, given 
k = 2,3. The author has further conjectured the following analogue of the Weyl 
inequality for hypergraphs, that is, if one can show that if two hypermatrices are 
“close” then they will have spectra which are “close”. Such a result would imply that 
the eigenvalues of K“) and (k — 1)!-!J* are asymptotically the same, because the 
difference of the entries of the aforementioned hypermatrices have norm converging 


to 0. 


Conjecture 8.11. (/11/) Suppose A and B are hypermatrices so that ||A — B|| <¢ 
for some norm ||-|| (or spectral radius) and é > 0. Then there is a bijection p between 
the eigenvalues (with multiplicity) of A and B and a function f with lim.) f(e) = 0 


so that |A — p(A)| < f(e) for each eigenvalue X. 


We present another approach to determining the eigenvalues of a random k-graph, 
by way of quasirandomness, which circumvents the need for a Weyl inequality. In- 
tuitively, a hypergraph is quasirandom if it has the same number of copies of any 
fixed subgraph as one would expect in a random hypergraph (where, in its simplest 
form, each edge is taken with probability 1/2). This idea was first introduced for 
graphs in [6] and was later extended to hypergraphs in [7]. In [7], Chung shows that 
a hypergraph is quasirandom if it has approximately the expected number of even 
partial octahedra (as described therein). One can restate this condition in terms of 
the coefficients, and perhaps the spectrum itself, by appealing to Theorem 3.27 to 
show that the linear combinations of subgraph counts appearing in the result are 
indeed forcing families for quasirandomness (see [26]). We can address this problem 


using the following definition. 
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Definition 8.12. Let 


bi(2) = (2? bn, (2))|zan/e = D0 che" 
n=0 


where D = deg(¢(H)). The characteristic power series of a random k-graph is 


(2) = dim Ejyjn[Oy,(z)]. 


The characteristic power series has a flavor of flag algebras (which should not come 
as a surprise as flag algebras are integral to the study of extremal combinatorics). We 
recommend [5] and [31] for readers interested in learning more about this topic. With 
this in mind, we adhere to the notation of flag algebras and denote P(H) simply as 


HA and further write G” to mean n disjoint copies of G. For example, 


Os(z) S140 4(-—P) 2" + (—263)2" (PF = 2C)2* 


+ (—2C5 + 2PoC3)z® + (—2C5 + 2C1P) + 4CZ — PR)2* +... 


is A 3 1 175 
ay a? BA GA 8a 


As an aside, we ask the following question: 
Question 8.13. Find a closed formula for 3. 


Famously, a forcing family of a graph is the edge and the 4-cycle so that, for @6(z), 
if c is ‘close’ to —1/2 and cj is ‘close’ to —1/4 then G is quasirandom. In particular, 
this implies (or more appropriately, forces) the coefficients of ¢% to be ‘close’ to the 
coefficients of 63. We have from [6] (and to some extent [26]) that the eigenvalues of 
a quasirandom hypergraph are asymptotically equal to those of the random graph so 
that the roots of ¢% are asymptotically equal to ¢3. The goal then is to apply this 
approach to ¢; for k > 3 (notice that the multi-set spectrum of the random k-graph 


is precisely the roots of ¢;). We can do so by answering the following two questions. 


Question 8.14. What are the forcing families of quasirandom hypergraphs? 
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More specifically, what is a small forcing family for quasirandom 3-graphs? An 
analogue of the graph case (i.e., edges and 4-cycles) is desirable. Such a result is 
achievable using flag algebras, and further, ¢; could be computed explicitly using 
flagmatic software for low-values of k. Finally, we need to guarantee that the roots 


of ¢; behave as expected under this limits. 


Question 8.15. If {H,}%2, is a sequence of k-graphs so that $3, (z) > ¢;(z) then 


does their spectrum similarly converge? 


Recall that this question is true for the graph case as the limit object of the 
sequence is the uniform graphon. The authors feel that a similar result would hold 


for the hypergraphon but the details will need to be verified. 
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